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1 .  Engineering  Significance  of  Phase  I  Study 

In  comparison  with  other  existing  fatigue  crack  growth  models,  the 
present  model  offers  the  following  features  of  engineering  significance. 

a)  The  mean  fatigue  crack  growth  model  is  developed  based  on  fundamental 
results  provided  by  physical  theory.  The  parameters  of  the  model  thus 
have  physical  significance  and  they  provide  new  insight  into  the 
crack  growth  process. 

b)  The  model  parameters  are  directly  related  to  the  simple  harmonic 
tensile  load  and  the  temperature.  Hence,  test  data  under  varying 
loading  conditions  and  varying  temperatures  can  be  used  for  parameter 
estimation.  Furthermore,  validation  of  the  model  using  existing  test 
data  shows  that  the  crack  growth  behavior  can  be  reliably  predicted 
based  on  two  or  at  most  three  suitably  chosen  tests.  Considerable 
savings  in  cost  and  time  associated  with  testing,  data  collection  and 
data  analysis  can  thus  be  realized. 

c)  The  model  indicates  that  testing  time  can  be  shortened  by  simply 
increasing  the  temperature  with  the  same  loads.  This  procedure 
eliminates  the  load  interaction  problem  when  using  the  standard  method 
of  accelerating  fatigue  and  FCG  tests  by  increasing  the  loads. 

d)  Inherent  variabilities  in  material  behavior,  load,  temperature,  and 
environment  can  be  incorporated  into  the  mean  model. 

From  an  engineering  viewpoint,  these  features  represent  an  important 
improvement  over  existing  and  more  phenomenologically  based  models.  It 
is  proposed  that  further  development  of  this  model  and  its  imp! ementational 
issues  can  be  addressed  in  Phase  II. 
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2.  Summary  of  Work  Done  in  Phase  I  Study 

The  following  tasks  were  specified  in  the  Phase  I  Proposal  and  in 
Contract  F49620-90-C-0059  funded  by  AFOSR: 


Task  i) 


Task  ii) 


Task  iii) 


Task  iv) 


A  deterministic  macro  model  of  damage  accumulation  for  fatigue 
crack  growth  (FCG)  will  be  constructed  from  the  micro  theory. 
Identify  the  parameters  in  the  macro  model  with  physical  and 
material  properties  of  engineering  significance  to  determine 
which  must  be  estimated  by  testing  and  which  can  be  estimated 
from  known  material  properties. 

Four  sets  of  existing  data  on  FCG  will  be  used  to  explore  the 
parameter  estimation  problem  in  order  to  determine  the  minimum 
amount  of  testing  needed  to  characterize  a  component  for 
prediction  purposes. 

Determine  how  variability  is  to  be  included  in  the  deterministic 
engineering  model  so  that  the  inherent  probabilistic  distributions 
present  in  the  phenomenon  can  be  included. 


These  tasks  and  more  were  completed. 
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Let  us  summarize  briefly  how  these  tasks  were  accomplished,  and 
what  was  learned  in  the  process.  Details  are  provided  in  the  references 
1  and  2. 

Figure  1  shows  the  crack  plane  on  one-half  a  thin  center  crack 
panel  under  axial  tensile  load;  aQ  is  the  initial  one-half  crack  length, 
a  is  the  current  one-half  crack  length,  and  £  is  a  generic  point  in 
the  intake  portion  of  the  one-half  crack  plane.  Let  us  introduce  some 
notation. 


aj  -  1/2  the  crack  length  considered,  j  =  1  ,  ...,  n, 
yj  =  aj/y  =  1/2  the  panel  width, 

Aaj  =  aj  '  aj-l'  ayj  •  Aaj' V 
w  =  panel  thickness, 

so  =  ^ar  ax->al  tensile  stress, 

=  sQ  +  ac  h(  wt)  >  0 
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h(x)  =  periodic  function  on  (0,2  tt)  ,  | h  (t  )|  _<  1 ,  mean  h  (t  )  =  0 

T  =  absolute  temperature  (Kelvin), 
t.  =  mean  time  to  reach  a..  At,  =  t-  -  t.  , . 

J  J  •*" 

2  (yj+yj-i)/2- 

Task  i):  A  deterministic  macro  model  of  damage  accumulation  for  fatigue 
crack  growth  (FCG)  will  be  constructed  from  micro  (reaction  rate)  theory. 

A  deterministic  macro  model  of  FCG  in  a  thin  center  crack  panel  of 
polycrystal ine  ductile  metal  under  an  axial  tensile  load  with  a  periodic 
time  variation  is  produced  from  a  micro  result  provided  by  reaction  rate 
theory.  Specifically,  the  macro  model  gives 

Atj  =  f(Ay.,  y j,  e,  z;  sQ,  As,  h(x),  T;  AG+,  Bi  or  \./\>  , 

where  g.  =  \,/vkT  is  the  growth  control  parameter,  k  =  Boltzmann's  constant, 

J  J 

h  =  Planck's  constant,  and  v  =  mean  number  of  bonds  per  unit  area.  Only 
one  damage  accumulation  observable,  namely  y,  is  present,  although  other 
possibilities  exist,  e  and  z  are  geometrical  constants  that  refer  to 
the  normal  stress  in  distribution  s(£)  on  (a  <_  £  £  lQ)’>  and  T  is  assumed 
constant.  We  note  that  to  predict  At- 

sJ 

the  geometrical  quantities  Ay,,  y,,  e,  z 

J  J 

must  be  selected, 

the  environmental  factors  5  ,  As  ,  h(t),  T 

must  be  obtained  from  the  operating  conditions,  and 

the  material  parameters  AG+,  6-  or  A,/v 

J  3 

must  be  estimated . 

The  model  is  new.  It  is  based  upon  physical  theory,  and  it  explicitly 
contains  qualities  of  engineering  interest.  Details  of  the  derivation 


are  given  in  Sections  2,  3,  and  4  in  Reference  [1]  and  in  Section  2  in 
[2].  Thus  Task  i)  is  completed  for  one  specimen  geometry  and  load. 
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Task  ii):  Identify  the  parameters  in  the  macro  model  with  physical  and 
material  properties  of  engineering  significance  to  determine  which  must 
be  estimated  by  testing  and  which  can  be  estimated  from  known  material 
properties. 

The  macro  deterministic  model  contains  explicitly  the  physical 
(environmental)  parameters  iQ,  as,  h(t),  and  T  that  are  of  direct  engineer¬ 
ing  significance,  i  ,  as,  and  h(t)  describe  the  load  considered  here,  and 
are  known  once  the  load  is  specified.  T  is  the  temperature  of  the  component 
or  specimen  and  is  known  once  the  environment  is  specified.  We  usually 
assume  T  is  constant,  although  this  is  not  necessary;  we  also  assume 
there  are  no  chemical  agents  in  the  atmosphere.  It  turns  out  from  our 
studies  that  h(t)  and  T  are  more  important  in  the  estimation  of  At.  than 

J 

previously  thought. 

The  macro  model  also  contains  the  material  parameters  aG+  and  g.  • 

vJ 

AG+  is  a  Gibbs  free  energy;  it  is  a  thermodynamic  property,  and  it  only 
depends  on  the  material  composition;  it  can  be  estimated  for  metals  from 
known  information  on  the  heat  of  sublimation  of  the  metal ic  elements  and 
the  composition.  ,3.  not  only  depends  on  the  material  composition  but 
it  also  depends  on  how  the  material  is  formed  (cast,  forged,  annealed,  etc.) 
into  components.  3.  is  correlated  with  microscopic  material  properties. 

6.  must  be  estimate  from  test  data  at  present.  Details  are  given  in 
Sections  2,  3,  and  4  in  [1],  and  in  Section  2  in  [2].  Task  ii)  is  thus 
accompl i shed. 
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Task  iii):  Four  sets  of  existing  data  on  FCG  will  be  used  to  explore  the 
parameter  estimation  problem  in  order  to  determine  the  minimum  amount  of 
testing  needed  to  characterize  a  component  for  predictive  purposes. 

Actually,  six  sets  of  FCG  data  were  used  to  explore  the  parameter 
estimation  problem.  The  parameter  we  estimate  is  where  s^  =  sQ+  as. 

In  each  set,  it  was  found  that  Bjig  could  be  expressed  in  terms  of  a 
linear  expression  in  8Q»  c  and  y:  =  BqS^  1-cy  • ) .  ^his  means  there 

are  three  material  parameters  aG+,  bo>  c  of  which  8Q  and  c  must  be  estimated 
from  test  data.  It  was  found  that  80  depends  on  the  load  whereas  c  does 
not.  bq  depends  on  T  in  a  known  manner  but  c  does  not  seem  to.  Tests 
at  three  different  loads  are  the  maximum  required  in  order  to  be  able  to 
estimate  B0  as  a  function  of  the  load.  Then,  it  is  possible  to  predict 
the  mean  time  t.  to  reach  the  crack  length  a.  for  a  range  of  loads. 

J  J 

We  also  found  that  it  is  possible  to  accelerate  tests  by  increasing 
T.  Thus,  we  have  found  the  minimum  amount  of  testing  required  to  character¬ 
ize  a  component  for  predictive  purposes. 

For  details,  see  Sections  3,  4,  and  7  in  [2]  and  Sections  5  and  6 
in  [1],  We  thus  have  completed  Task  iii). 

Task  iv):  Determine  how  variability  is  to  be  included  in  the  deterministic 
engineering  model  so  that  the  inherent  probabilistic  distributions  present 
in  the  phenomenon  can  be  included. 

We  have  found  that  variability  can  be  introduced  at  the  macro  level 
in  the  deterministic  model  by  assuming  that  the  work  done  in  bond  breaking 
on  the  Ay,  interval,  8. s',  has  a  random  factor  Q..  The  data  sets  are 
employed  to  estimate  the  moments  and  distributions  of  the  Q,.  What  is 

<J 

remarkable  is  that  the  distributions  of  the  standardized  Q.  are  very 

J 
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close  to  one  another  and  do  not  seem  to  depend  on  the  load  or  interval 
length.  We  have  no  explanation  for  this  remarkable  fact  based  on  this 
macro  method  for  introducing  variability.  However,  this  approach  to 
variability  does  not  enable  us  to  determine  if  the  distributions  of  the 
Q.  depend  on  T.  The  tests  needed  for  parameter  estimation  in  the 
deterministic  model  suffice  to  estimate  the  parameters  in  the  probabil¬ 
istic  model  for  predictive  purposes  providing  T  does  not  change. 

See  Sections  5  and  6  in  [2]  for  details.  Task  iv)  is  thus  completed. 
Persual  of  [1]  and  [2]  makes  clear  that  more  must  be . done  to  validate 
the  current  model,  to  further  develop  the  current  model  so  that,  for 
example,  it  applies  to  different  crack  geometries  and  loads,  and  to  remove 
some  of  the  restrictive  assumptions  employed  in  the  current  model. 
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CUMULATIVE  DAMAGE  MODEL  FOR  MEAN 
FATIGUE  CRACK  GROWTH  BASED  ON  THE 
KINETIC  THEORY  OF  THERMALLY 
ACTIVATED  FRACTURE 

F.  KOZINt  and  J.  L.  BOGDANOFFJ 

tPolytechmc  University  Route  1 10.  Farmingdaie,  NY  (1735.  U  S. A. 

^School  of  Aeronautics  and  Astronautics,  Purdue  University,  W.  Lafayette.  IN  47907,  U  S  A. 

Abstract — The  objective  of  this  paper  is  to  show  how  the  physical  concepts  due  to  Tobolsky-Eyring 
for  the  description  of  atomic  (micro)  bond  fracture  can  be  extended  to  the  macro  problem  of  mean 
fatigue  crack  growth  which  is  of  engineering  interest.  The  parameters  of  the  model  have  physical 
meaning.  Parameter  estimation  is  studied  with  the  aid  of  four  data  sets  on  fatigue  crack  growth. 
Because  of  the  fundamental  physical  concepts  employed,  the  model  provides  new  insights  into  the 
fatigue  crack  growth  process. 


1.  INTRODUCTION 

The  purpose  of  this  paper  is  to  present  an  initial  study  of  a  mean  fatigue  crack  growth  (FCG)  model 
based  upon  the  kinetic  theory  of  thermally  activated  atomic  bond  fracture  due  to  Tobolsky-Eyring[l  ]. 

We  start  from  a  basic  expression  for  the  mean  rate  of  single  atomic  bond  breaking  that  has 
its  origin  in  the  statistical  theory  of  thermally  activated  bond  fracture.  A  brief  discussion  of  the 
applicability  of  the  expression  is  given. 

From  the  basic  expression,  which  applies  at  the  micro  level,  we  derive  a  differential  equation 
for  mean  crack  growth  in  a  center  crack  panel  or  WOL  specimen  under  tensile  load  that  applies 
at  the  macro  level.  This  is  a  major  contribution  of  the  paper. 

Employing  certain  simplifying  assumptions,  we  then  obtain  a  formula  for  mean  crack  growth 
whose  parameters  have  physical  meaning.  Some  of  these  parameters  are  directly  related  to  the  load; 
others  refer  to  material  properties.  One  parameter  can  frequently  be  estimated  from  tabulated 
physical  data;  however,  the  other  unknown  parameters  must  be  estimated  from  experimental  data. 

The  statistical  analysis  available  for  four  sets  of  high  replication  data  provide  the  basic 
information  on  mean  crack  growth  as  a  function  of  time  with  which  the  remaining  parameters  can 
be  estimated.  These  estimates  provide  new  insights  into  the  crack  growth  process  under  simple 
harmonic  tensile  loading,  including  the  effect  of  load  shedding.  This  is  a  second  major  contribution. 

A  third  major  contribution  is  the  fact  that  the  mean  model  presented  can  predict  mean  fatigue 
crack  growth  behavior  in  a  component  under  arbitrarily  varying  simple  harmonic  tensile  loading 
at  any  temperature  based  on  two  or  at  most  three  suitably  chosen  tests. 

Finally,  we  discuss  the  results,  and  offer  suggestions  for  future  study  part  of  which  will  be 
presented  in  a  subsequent  paper. 

We  wish  to  emphasize  at  this  point  that  the  model  we  develop  is  different  from  other  mean 
models.  We  are  concerned  with  cumulative  damage  in  which  crack  length  is  a  macro  observable 
and  a  measure  of  macro  damage  accumulation,  starting  from  a  basic  expression  for  atomic  bond 
breaking  which  applies  at  the  micro  level.  Earlier  models  do  not  have  this  focus.  Further,  this  model 
has  no  contact  with  the  concepts  usually  employed  in  fracture  mechanics. 

2.  BOND  BREAKING  WITH  A  SINGLE  ENERGY  BARRIER  AND  NO 

BOND  MENDING 

The  mean  rate,  k„.  of  single  atomic  bond  breaking,  when  there  is  only  a  single  energy  barrier, 
no  bond  mending,  and  a  tensile  external  force  present,  is,  according  to  physical  theory  as  developed 
by  Tobolosky-Eyring[l]  (see  also  p.  338  in  ref.  [2]),  given  by 
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where 


k  =  Boltzmann’s  constant, 

T  =  absolute  temperature, 
h  -  Planck’s  constant, 

AG  *  =  free  energy  of  activation, 

A  =  average  distance  over  which  force  acts  in  bond  breaking, 

/  =  average  force  acting  on  a  single  bond  during  bond  breaking. 

Let  us  briefly  discuss  the  meanings  of  the  last  three  quantities,  and  the  form  of  the  exponent  in  eq.  ( 1 ). 

AG*  and  A  are  material  parameters.  Tobolsky-Eyring[l,  2]  interpret  the  free  energy  AG*  as 
the  height  of  the  single  energy  barrier  that  must  be  overcome  in  order  to  break  a  bond.  AG  *  also 
may  be  thought  of  (Zhurkov[3])  as  the  binding  energy  of  atoms  in  metals  and  the  energy  of 
breakage  of  the  chemical  bonds  in  chain  macro  molecules  in  polymers. 

The  atomic  bond  force  / depends  on  the  external  tensile  force  acting.  We  shall  assume  below 
that  /  may  be  approximated  by  a  term  proportional  to  the  stress,  as  was  done  in  refs  [1—3].  Thus, 
eq.  (1)  says,  as  noted  in  ref.  [4],  that  “the  free  energy  of  activation  is,  in  effect,  reduced  linearly 
(in  stress)  in  a  manner  identical  to  that  proposed  by  Eyring[5]  for  molecular  interpretation  of  liquid 
viscosity”.  Further,  the  extensive  experimental  studies  of  Zhurkov[3]  also  demonstrate  that  a  linear 
dependence  on  stress  in  the  exponent  of  eq.  (1)  is  required  to  describe  his  data  on  time  to  failure 
under  constant  stress. 

The  second  material  parameter  is  A.  A  is  the  average  distance  over  which  the  force / acts  during 
bond  breaking.  In  ref.  [4],  it  is  regarded  as  the  mean  separation  between  the  equilibrium  positions 
of  minimum  force  potential.  It  is  clear  from  the  form  of  the  exponent  in  eq.  (1)  that  A  is  inversely 
proportional  to  the  strength  of  the  material.  A  will  be  influenced  by  how  the  material  is  formed 
into  components,  and  by  work  hardening  or  softening.  Zhurkov[3]  (replace  his  y  by  A)  gives  some 
very  interesting  observations  on  the  meaning  of  this  parameter.  For  example,  a  direct  correlation 
has  been  found  between  A  and  the  disorientation  of  the  structure;  the  more  disorientated  the 
structure  in  a  metal,  i.e.  the  higher  the  dislocation  density  at  the  slip  planes,  the  smaller  the  strength 
coefficient  A,  and  hence  the  higher  the  metal  strength,  etc.  Thus,  A  is  an  important  parameter. 

In  all  events,  the  form  of  eq.  (1)  has  substantial  support  in  experiments  on  a  variety  of  materials 
including  polycrystalline  metals  and  polymers. 

Equation  ( 1 )  applies  at  the  micro  scale;  i.e.  it  gives  the  mean  rate  of  single  bond  breaking  under 
the  conditions  stated.  Our  concern  with  failure  (fracture)  is  at  the  macro  level.  Thus,  it  is  necessary 
to  derive  from  eq.  (1)  a  macro  model  of  failure.  It  is  worth  noting  at  this  point  that  a  comparison 
of  several  macro  rate  reaction  models  for  fracture  in  solids  is  given  in  ref.  [4];  this  comparison 
indicates  that  the  models  based  on  bond  breaking  rather  than  on  bond  slipping[6]  are  consistent 
with  experimental  results  obtained  on  fracture  in  polymers  and  polycrystalline  metals.  However, 
in  these  models  of  fracture  the  growth  of  a  single  fatigue  crack  is  not  considered,  which  is  of 
engineering  interest. 

Our  approach  will  be  different  from  those  in  refs  [1-3]  in  that  the  growth  of  a  single  crack 
is  of  major  importance,  i.e.  we  are  concerned  with  a  cumulative  damage  process  in  which  there 
is  an  observable. 


3.  DERIVATION  OF  A  DAMAGE  ACCUMULATION  MODEL 

We  take  the  simplest  approach  in  this  initial  study.  To  fix  ideas,  let  us  consider  a  thin 
center-crack  panel  under  axial  tension.  Let  iv  denote  the  thickness  of  the  panel,  /„  the  half-width, 
and  F(t )  the  axial  tensile  force. 

Consider  the  cross-section  of  the  pane!  formed  by  a  plane  perpendicular  to  the  axis  of  the  panel 
containing  the  crack  which  we  assume  lies  in  this  cross-section.  Let  a  denote  the  half  mean  crack 
length. +  The  panel  is  assumed  thin  so  that  we  can  assume  the  crack  tip  is  straight  and  perpendicular 
to  the  panel  face.  Let  ^  denote  the  distance  along  the  axis  of  the  cross-section  from  the  centerline 


rWhen  we  speak  of  crack  length  we  always  mean  "mean  "  crack  length  from  a  distribution  at  lime  r  even  when  this  is 
not  explicitly  stated 
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to  a  generic  point  in  the  portion  of  the  cross-section  ahead  of  the  crack  tip,  i.e.  a  <  c,  /0  (see 
Fig.  1).  We  observe  that  WOL  geometry  is  the  same  as  one-half  a  center  crack  panel,  although 
stress  distributions  are  different. 

Let  v  denote  the  average  number  of  atomic/nuclear  bonds  per  unit  cross-section  area.  Then 

vw  d£  =  average  number  of  such  bonds  in  the 

interval  (s,  <;  +  d£)  in  the  unbroken  portion. 

If  we  assume  that  no  appreciable  bond  breaking  occurs  ahead  of  the  crack  tip,  it  follows  that 

% 

2vw  d£  =  total  average  number  of  unbroken  bonds 
Ja  in  the  unbroken  portion  of  the  cross-section  (2) 

=  2  vw(lo-a). 

This  is  an  important  step  because  it  connects  the  macro  variable  crack  length  a  with  the  micro  bond 
breaking  process.  Since  a  is  the  only  quantity  that  changes  with  time,  eq.  (2)  implies  that  bond 
breaking  is  concentrated  in  the  region  of  the  crack  tip. 

Let  s(t,  £)  denote  the  normal  stress  at  the  point  £  in  the  unbroken  cross-section  at  time  t.  We 
assume  s(t ,  =  s(t,  —c,)  and  that  this  stress  does  not  vary  across  the  cross-section.  We  note  that 
s(t,  £)  is  an  average  macro  concept[18]. 

The  normal  force  acting  on  the  interval  (£,  <;  +d{)  is 

*-$(/,£)  d?. 


If  we  divide  this  average  force  by  the  mean  numbers  of  bonds,  wv  d£,  in  this  interval,  we  obtain 
an  approximation  for  the  force  /  acting  on  a  bond  at 


(3) 

We  now  replace  eq.  (1)  by 

V 

»  ejp  AG*  -  AJ(/.  4)/V 

Kb(f,0  =  ye  kT  , 

(4) 

or,  with 

kT  ~  .  A 

a  =  Te  • 

(5) 

Kb(t,f)  =  ae*<'{\ 

(6) 

which  is  the  mean  rate  of  single  bond  breaking  at  point  c  at  time  t. 


Fig  1.  Geometry  of  one-half  of  a  center  crack  panel  with  some  notation. 
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The  macro  average  rate  of  bond  breaking  on  the  unbroken  portion  of  the  cross-section  is  thus 


% 

J“ 


2  wvKbd<;=2wv  %  e*1'-4’  di;. 


(7) 


Let  us  set  this  equal  to  the  negative  derivative  with  respect  to  time  of  the  total  average  number 
of  bonds,  2 vw(l0  —  a),  in  the  unbroken  portion  of  the  cross-section,  obtaining 


da 

dt 


i  e,,l'-{>  d<;. 


(8) 


Neglecting  dynamical  effects,  we  have  from  equilibrium  at  time  r. 


F{t)  =  2w 


s(t,  c)df. 


(9) 


We  note  that  in  eqs  (8)  and  (9)  s(t,  <;)  occurs  under  an  integral  sign  with  respect  to  This  means 
that  the  right  hand  sides  of  eqs  (8)  and  (9)  will  not  be  overly  sensitive  to  how  s(t,  £)  depends  on 
since  s(t,q)  is  a  bounded  function.  Further,  we  note  that  eq.  (8)  is  not  of  the  form  usually 
encountered  in  fracture  mechanics;  it  is  new  and  has  been  derived  from  physical  theory.  Since  eqs 
(8)  and  (9)  are  basic  in  what  follows,  it  is  important  at  this  point  to  make  some  comments  on  their 
derivation. 

Simply  put,  to  produce  eq.  (8)  we  evaluate  the  mean  number  of  unbroken  bonds  at  time  t  based 
only  on  geometry  (see  eq.  2)  and  equate  the  negative  derivative  with  respect  to  time  of  this  number 
to  an  expression  for  the  mean  rate  at  which  bonds  are  being  broken,  as  determined  by  the  kinetic 
theory  of  thermally  activated  fracture  (see  eq.  7).  The  time  dependent  quantity  in  eq.  (2)  is  the  crack 
length  a{t)\  thus  damage  accumulation  in  this  expression  is  measured  by  the  macro  quantity  a(t), 
with  rate  of  accumulation  proportional  to  a(t).  On  the  other  hand,  the  rate  of  damage 
accumulation  according  to  eq.  (7)  depends  on  s(t,  <J)  as  well  as  on  a(/).  s(t,  £)  not  only  depends 
on  a(t)  but  also  on  other  quantities  such  as  plastic  zone  size,  for  example.  We  note  that  eq.  (8) 
can  allow  the  possibility  that  additional  macro  observables  may  be  present. 

As  stated,  eq.  (9)  is  an  equilibrium  result.  This  means  that  we  assume  the  frequency  content 
of  F(t)  is  well  below  the  lowest  natural  frequency  of  the  specimen  that  can  be  excited  by  axial 
loading.  However,  if  dynamical  results  must  be  included,  5(r,  ?)  will  be  related  to  F(t)  by  a  more 
complex  relation  than  that  given  by  eq.  (9). 

We  have  approximated /(see  eq.  3)  by  means  of  the  average  macro  quantity  s(t,  <;)■  Thus,  we 
are  employing  a  smoothed  version  of /  However,  in  a  real  material,/ as  well  as  A  will  change  from 
bond  to  bond  due  to  variability  in  material  properties  as  noted  by  Zhurkov[3].  It  follows  that  when 
such  variability  is  to  be  taken  into  account,  we  must  re-examine  our  evaluation  of  A /. 

x  and  P  contain  the  temperature  T.  In  the  constant  stress  tests  of  Zhurkov[3],  the  temperature 
was  held  constant.  However,  T  may  increase  when  F(t)  is  cyclic,  as  is  well  known  when  the 
frequency  is  high.  Hence,  T  may  be  time  dependent. 

Next,  p  contains  the  material  parameter  A.  Zhurkov[3]  as  well  as  others[l,2]  assumed  A  was 
a  constant  in  time  to  failure  under  static  tests,  and  this  was  not  contradicted  by  the  experimental 
results.  However,  where  damage  accumulation  is  being  studied  under  cyclic  load,  this  may  not  be 
the  case. 

Finally,  we  observe  that  eqs  (8)  and  (9)  apply  to  a  particular  geometry,  i.e.  to  either  a  center 
crack  panel  or  to  a  WOL  specimen.  Other  geometries  of  crack  development  will  require  a  different 
formulation.  Such  possibilities  will  not  be  considered  in  this  initial  study. 


4.  SOME  SIMPLE  MACRO-MODELS  OF  DAMAGE  ACCUMULATION 

At  this  initial  stage  in  our  investigation  of  what  can  be  learned  from  eqs  (8)  and  (9),  it  is 
appropriate  to  confine  attention  to  simple  situations.  To  that  end.  let  us  first  assume 


J(f?)  =  i«). 


(10) 
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which  means  that  the  stress  is  uniform  over  the  unbroken  portion  of  the  panel  cross-section.  We 
first  find  from  eq.  (9) 


i(f)  = 


F(t) 

2  w(l0-a)' 


(ID 


From  eq.  (8),  we  then  obtain 


da 

dt 


2 


f  V'^ds  =  2(/0-a)e*,,) 


(12) 


or 


tm 


The  substitutions 


da 

—  =  a  (4  —  a)  e2“',,°-a). 
dr 


y=T’  77  =  so(e), 
‘o  2  w/0 


where  50(r)  is  the  far  field  axial  stress,  change  eq.  (13)  to 

dv 

^=«(1  -y)e‘~\ 


(13) 


(14) 


(15) 


Different  forms  of  this  equation  appear  in  refs  [2, 4],  where  it  is  assumed  s0(r)  =  constant,  i.e.  the 
constant  load  case.  When  the  load  or  stress  is  constant,  eq.  (15)  can  be  integrated  to  yield  a  time 
to  failure  expression(2],  and  it  is  this  expression  that  not  only  provides  an  excellent  description  of 
the  time  to  failure  data  due  to  Zhurkov(3]  at  different  constant  loads  and  different  temperatures, 
but  it  also  provides  an  accurate  predictor  of  time  to  failure  where  the  load  and/or  temperature  are 
changed  from  the  values  at  which  a  and  /?  are  initially  estimated. 

Zhurkov  actually  does  not  start  from  eq.  (15).  Rather  he  employs  a  time  to  failure  expression 
obtained  by  fitting  an  exponential  formula  to  the  extensive  data  that  he  collected.  The  empirical 
expression  he  obtained  can  be  obtained  from  the  integrated  form  of  eq.  (15)  by  using  suitable 
approximations.  For  a  discussion  of  this  point,  see  ref.  [5]. 

In  refs  [6,  7],  s(t)  is  assumed  simple  harmonic  about  a  mean  value  (>0).  Equation  (8)  is  then 
numerically  integrated,  and  the  results  compare  well,  in  the  example  considered,  with  the  single 
specimen  test  result.  In  ref.  [7],  they  also  consider  this  case  employing  a  different  method. 

Equation  (15)  is  the  simplest  form  we  can  obtain  from  eqs  (8)  and  (9).  It  has  one  observable 
or  measure  of  damage  accumulation,  namely  y.  However,  an  objection  to  eq.  (15)  can  be  raised. 
It  stems  from  the  fact  that  the  stress  is  assumed  uniform,  i.e.  independent  of  c,  which  means  there 
is  no  physical  reason  why  there  cannot  be  significant  bond  breaking  well  ahead  of  the  crack  tip 
although  we  assumed  in  eq.  (2)  the  contrary.  To  partially  overcome  this  objection  let  us  next 
assume. 

s(?,O  =  *(Og(0,  (16) 

where  the  stress  factor  j(<J)  only  depends  on  £  and  g(t)  is  a  multiplicative  time  factor.  In  view  of 
the  fact  that  the  normal  stress  is  high  in  the  irreversible  bond  breaking  region  and  then  decreases  to 
an  elastic  stress  as  we  move  away  from  the  crack  tip,  a  more  realistic  assumption  for  ^(^)  is  either 

*(<?)  =  <V  a<f<a  +  A,  (17a) 

=  <re(  <  <7„),  a  +  A  <  £  $  /„, 


or 

j(<;)  =  <*p  +  — r—  Act  -  a  <!;  <  a  +  A,  (17b) 

A  A 

=  ctc,  a  +  A  <  c  ^  /„, 

where  A  <  /„  -  a  and  0  <  Aa  <  <7„  -  ac.  Equation  (17a)  gives  two  constant  stress  regions  as  does 
the  second  region  in  eq.  (17b).  However,  in  the  region  a  ^  ^  a  +  A,  eq.  (17b)  permits  the  stress 
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to  vary  linearly  from  crp  +  A<r  to  cp  —  A<r  >  <jc  with  mean  value  ffp.  Thus,  in  eq.  (17b),  the  stress 
at  the  crack  tip  is  the  largest  in  the  unbroken  cross-section,  which  means  that  the  maximum  rate 
of  bond  breaking  occurs  at  the  crack  tip  even  if  Aa  is  small.  If  Ac  =  0,  eqs  (17a)  and  (17b)  are 
the  same. 

We  find  from  eq.  (9),  using  either  eq.  (17a)  or  (17b),  that  the  dimensionless  time  factor  is 


gU)  = 


Sp(t)l<Tc 

(1  -y -&  +r„S)' 


(18) 


where  S  =  A/70,  a p/{Tc  =  ra=  1  +  £,  and,  as  before,  s0(t)  =  F(t)/2wl0  is  the  far  field  stress.  The 
substitution  of  eqs  (18)  and  (17a)  into  eq.  (8)  produces 

j  ft,  jqI')  Mo 

—  =  x[d  e1  - y -**'„>  q-(l  —  j.-  —  <5)e‘  (19a) 

at 


if  eq.  (17b)  is  employed,  eq.  (19a)  is  replaced  by 


dy 

d7 


=  otS 


+') 

gl  -  y  -  6  0q6 


sinhl 


Ps0(t)Ao'  \ 
1  -  y  -  3  +  r„3  J 
/3s0(t)Acr'  \ 


,1  -y-S  +  r,6  j 


food 

+  a(l  -y  —  3  +  r„<5)e|->'~'i  +  r«‘>, 


(19b) 


where  Act'  =  Aa/ac.  If  S  =0,  both  reduce  to  (15). 

Equation  (19b)  is  approximately  equal  to  eq.  (19a)  for  A  o’  1.  Hence,  we  shall  use  eq.  (19a) 
in  subsequent  work  since  it  is  the  simpler  of  the  two  from  the  parameter  estimation  point  of  view, 
and  we  do  not  have  to  consider  A  a'. 

We  note  that  in  eq.  (19a)  there  are  three  observables  or  measures  of  damage  accumulation, 
namely,  y,  <5  and  r„  or  e.  We  must  either  obtain  expressions  for  the  joint  evolution  of  y,  3  and  £. 
or  we  must  assume  that  3  and  e  depend  on  y  in  some  explicit  manner  if  y  is  the  only  one  observable. 
We  shall  see,  after  confronting  data,  that  for  these  data  the  major  portion  of  atomic  bond  breaking 
does  not  occur  in  the  region  a+A<(;^/0if£>  0.20. 


5.  APPLICATION  OF  EQS  (15)  AND  (19)  EXPERIMENTAL  DATA  ON  MEAN 

CRACK  GROWTH 


The  purpose  of  this  section  is  to  confront  the  model  defined  by  eqs  (15)  and  (19)  with  data. 
We  shall  consider  four  sets  of  data  on  FCG  under  simple  harmonically  varying  tensile  load: 

(a)  ASTM  E24.04.01  round  robin  data,  10Ni-8Co-lMo  steel,  24  replications(8], 

(b)  Virkler,  Hillberry,  Goel  (VHG)  data,  2024-T3  aluminum,  68  replications^,  10]. 

(c)  Ichikawa,  Hamaguchi,  Nakamura  (IHN)  data,  2024-T3  aluminum,  30  replications[l  l]. 

(d)  Ghonem,  Dore  (GD)  data,  7075-T6  aluminum,  60  replications[10,  12], 

Equations  (15)  and  (19)  refer  to  mean  crack  length  a  or  its  dimensionless  version  y.  For 
parameter  estimation,  we  require  an  estimate  of  how  a  or  y  varies  with  cycle  number  or  time.  Mean 
crack  length  is  obtained  by  averaging  over  the  a-values  of  the  sample  functions  (of  a  vs  time)  at 
a  fixed  time.  Available  data  analysis[10]  is  for  mean  time  to  reach  a  given  crack  length  that  is 
obtained  by  averaging  over  the  time-values  of  the  sample  functions  at  a  fixed  crack  length.  In  either 
case  a  y  vs  t  curve  is  obtained  in  which  the  term  "mean”  is  applied  to  only  y  or  to  t.  Examination 
of  the  VHG  data  presented  in  ref.  [10]  reveals  that  the  two  curves  are  close  to  one  another  (see 
Appendix).  Thus,  for  the  y  vs  t  curve,  we  shall  use  the  mean  time  vs  y  curve  which  is  sufficiently 
close  to  the  mean  y  vs  t  curve  for  our  purposes. 

Equations  (15)  and  (19)  contain  the  two  physical  parameters 


i  = 


fi 


/. 

v. kT " 


(20) 
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AG  *  can  be  estimated  in  one  of  two  ways.  Zhurkov  has  estimated  A G  *  from  simple  tests  on  time 
to  failure  under  static  load,  however,  Zhurkov[3]  has  also  shown  that  AG  *  can  be  obtained  from 
the  heat  of  sublimation  in  metals[  1 3]:  thus,  there  are  two  methods  for  estimating  AG  and,  hence 
for  estimating  a.  cannot  be  estimated  from  existing  physical  data  since  it  depends  on  the  test 
conditions  plus  the  manufacturing  process.  It  follows  that  /?  must  be  estimated  from  test  data.  We 
note  that  Zhurkov’s  simple  tests  for  estimating  /?  under  static  load  may  give  values  that  do  not 
apply  with  cyclic  loads. 

The  data  in  the  four  cases  consist  of  a  vs  n  (cycles)  from  which  we  obtain  the  values  of  A r, 
corresponding  to  Aa,  or  A y,.  We  shall  use  (A/,,  Ay,)  to  estimate  /J  in  eq.  (15)  for  the  first  data  set 
and  in  eq.  (19)  for  all  four  data  sets.  There  are  several  reasons  for  doing  this:  firstly,  it  is  important 
to  determine  whether  or  not  /?  depends  on  y  =a/lQ  or  time  /;  secondly,  if  /?  depends  on  these 
quantities,  it  is  important  to  determine  if  that  dependence  can  say  something  about  the  FCG 
process  and  its  modeling;  finally,  it  is  important  to  determine  what  must  be  known  in  order  to  make 
accurate  predictions. 


(a)  ASTM  £24.04.01  round  robin  data[\Qi\ 

The  material  is  10Ni-8Co-lMo  steel.  The  WOL-type  test  specimens  0.250  in.  thick  were  cut 
from  the  same  sheet  of  material.  The  load  employed  is 


F(t)  =  605  +  495  sin  cot  (lb.). 


(21) 


where  co  =  lOn/s.  Tests  were  conducted  at  different  laboratories.  From  the  results  obtained,  we 
selected  24  specimens  as  usable  for  our  statistical  aualysis[14].  A  load  shedding  program  was 
employed  to  produce  the  initial  crack  length,  which  differed  among  the  laboratories.  It  is 
appropriate,  therefore,  to  take  the  initial  crack  length,  a0  =  1 .0  in.  With  A  a  =  0.05  in.  and 
/0  =  2.55  in.  The  (/,  a)  data  are  available  at  20  points  for  purposes  of  parameter  estimation. 

The  material  is  steel.  If  we  regard  the  material  to  be  iron,  then  according  to  ref.  [15] 
AG  *  =  94,000  cal/mol.  We  take  T  =  300  K  as  a  representative  temperature,  and  find 

a  =  1.9655  x  10-5</s.  (22) 


The  parameter  /J,  must  be  estimated  from  the  (A/,,  Ay,)  data. 
Let  us  approximate  eq.  (15)  by  the  difference  form 


where 


£ih-  (|  4-  Ai'  sin  t) 

e1*^  dt 


Ay,  =  A a,//0  (y'th  y-increment), 

9i  =  (mid  value  for  y, » y, «. i ). 

Atj  -  times  increase  of  y}  by  A yr 
^J  =  vkf  °r  =  v^T  ’ 


(23) 


s0  =  constant  far  field  mean  stress. 

As  =  time  varying  far  field  stress  amplitude, 
A  s'  =  As  t sn ,  and 


r  =  cot. 


For  this  initial  study,  we  will  make  the  computational  simplifications  of  replacing  the  sine  wave 
form  by  a  rectangular  wave  form: 

-  for  0  <  t  <  7t, 
n 

_  2 
K 


for  7r  <  r  ^  2 n. 


(24) 
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Then 


where 


_t_ 

2k 


llJl.  (|  Ai'  «n  r) 

el-*i  dr 


(25) 


s'  =  £o  1  + 


(26) 


and  where  the  term  arising  from  the  second  part  of  eq.  (24)  is  neglected  since  in  the  data  analysis 
it  will  be  found  that  this  term  is  very  small  relative  to  the  contribution  from  the  first  part  of  eq.  (24). 
Equation  (23)  is  thus  replaced  by 


k  — 

Ay,=  2^  ->?;)A/>e|-  f/. 


(27) 


0T.\M£a. 


Since  Ay,  a,  y}  and  At,,  are  known,  we  can  estimate  0,io  from  this  equation. 

Figure  2  shows  the  estimates  as  a  function  of  y,  (or  d,)  and  When  the  estimates  are 
plotted  against  tJ  (Fig.  2a)  a  decreasing  nonlinear  function  is  obtained.  On  the  other  hand,  when 
they  are  plotted  against  y,  (Fig.  2b,  top  line),  the  straight  line  character  of  these  estimates  is 
remarkable.  Since  it  is  easier  to  work  with  a  straight  line  than  a  nonlinear  function,  we  prefer  the 
la The  least  square  linear  fit  to  these  estimates  is 

fiyS'o  =  116.921(1  -0.98788y),  (28) 


where  fiys'0  is  dimensionless.  Before  discussing  some  of  the  implications  of  this  result,  let  us  use 
the  form  in  eq.  (19a). 

Equation  (19a)  contains  the  two  additional  parameters  r„  and  5.  Let  us  replace  /•„  by  1  +  c. 
We  know  the  plastic  region  increases  in  length  as  the  crack  grows.  In  this  initial  study,  we  make 
the  simple  assumption  that  <5  =  zy,  0  <  z  <  1,  recognizing  that  this  is  an  assumption  of  importance 
and  it  will  be  discussed  below.  Thus,  we  introduce  the  constraint 


r,=  l+£,  <5=zy,  (29) 

where  c  and  z  constant.  We  follow  the  same  procedure  as  we  employed  in  obtaining  (27)  to  find 
from  (19a) 

S,<  I  -t-tlio  _  Si  Si> 

Ay,  =  -  zy,  At,  e1  _<l  +  -[!—(!+  z)y,]  At,  e1  (30) 


Fig.  2.  (a)  13, s’v  estimate  from  the  ASTM  data  plotted  against  time  (b)  0,s’o  estimates  plotted  against 
v,  based  on  eqs  (15)  and  (19).  Ay  =  0.05/2.55. 
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Fig.  3.  Comparison  of  estimated  and  observed  v  vs  t  points  for  ASTM  data,  where  estimated  curve 
employs  eq.  (27)  and  uses  /},s'0  from  eq.  (32)  rather  than  data. 

Let  us  take  t  =  1  and  z  =  0.25  as  representative  for  our  present  purposes.  We  quickly  find  on 
confronting  the  data  that  the  second  term  is  orders  of  magnitude  smaller  than  the  first  term  on 
the  right  hand  side  of  eq.  (30).  Hence,  eq.  (30)  becomes  in  practice 

_ 

Ay,  =  0.1 25a>:;  Ary  e1  ~  075  >>.  (31) 

Figure  2(b)  (see  lower  line)  also  shows  the  estimates  ^s'0  vs  y t.  The  straight  line  character 
of  these  estimates  still  obtains.  The  least  squares  estimate  of  this  line  obtained  from  the  estimates 
is 


pys'0  =  59.2759  (1  -  0.73625y).  (32) 

Thus,  regardless  of  whether  we  use  eqs  (15)  or  (19a)  the  s’Q  estimates  vs  y,  can  be  accurately 
approximated  by  a  straight  line.  What  does  this  result  say? 

Consider  ps'0  -  (ks'0/vkT)  ■  s'0  is  a  constant  as  are  v,  k  and  by  assumption  T.  Thus,  if  ps'0 
decreases  with  y,  this  means  that  k  decreases  with  y  also.  2  determines  the  strength  of  the  material, 
where,  as  discussed  in  connection  with  eq.  (1),  when  k  is  high  means  that  the  material  is  weaker 
than  when  k  is  low.  Figure  2  shows  that  k  is  decreasing  linearly  with  >■;  this  means  the  material 
(just  ahead  of  the  crack  tip)  is  getting  stronger  as  the  crack  length  increases;  this,  in  turn,  means 
that,  since  there  is  no  external  source  acting  to  produce  a  material  change,  the  material  is  work 
hardening  at  a  constant  rate  as  y  increases.  What  we  learn  from  this  model  is  that  the  rate  at  which 
work  hardening  is  taking  place  is  constant  as  a  function  of  y,  and  this  is  a  new  way  to  demonstrate 
that  fact.  The  same  results  hold  whether  we  use  eq.  (15)  or  (19a)  which  is  also  remarkable. 
Further  comments  will  be  delayed  until  remaining  data  sets  are  examined. 

Figure  3  plots  the  mean  crack  growth  curve  as  observed  and  one  predicted  from  (27)  with  s’0 
estimated  from  (32)  rather  than  from  the  data.  The  agreement  is  all  one  could  expect. 

(. b )  VHG  data[ 9,  10.  14] 

The  material  is  2024-T3  aluminum  alloy.  The  specimen  is  a  center  crack  panel  558.8  mm  long 
with  a  half-width  /„  =  76.2  mm  and  thickness  of  2.54  mm.  The  load  is 

F(t)  =  3150  -I-  21 10  sin  cut  (lb.)  (33) 

with  w  =  40  k  rad,  s,  68  replications  were  used.  Data  were  recorded  as  cycles  to  reach  a  given  crack 
length.  A  load  shedding  program  was  employed  to  produce  the  initial  crack  of  an  =  9  mm;  the 
maximum  crack  length  employed  was  49.8  mm.  Three  crack  length  intervals  A  a  were  employed, 
namely,  0.2  mm.  0.4  mm  and  0.8  mm,  /„  =  76.2  mm  and  (t,  a)  is  available  at  different  numbers  of 
points. 


♦Here  we  use  AC'  =  93.040  cal  mo!  which  lakes  into  account  the  fact  that  the  material  is  not  all  iron  This  gives 
i  =  1.94547  *  |0-'*  s 
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Fig.  4.  (a)  P,s'0  estimates  from  VHG  data  plotted  against  time  t .  (b)  PJ'0  estimates  plotted  against  >■,. 

A y,  =  1/76.2. 

The  data  analysis  in  refs  [10,  14]  estimated  the  mean  time  to  reach  a  specified  crack  length. 
As  mentioned  earlier,  we  shall  use  these  available  data  to  estimate  mean  crack  length  at  a  specified 
time  since  the  difference  is  small  between  the  two  cases. 

Zhurkov[3]  found  from  experimental  data  that  AC*  =  53,000 cal/mol  for  polycrystalline 
aluminum,  whereas  [14]  gives  55,000  cal/mol  for  aluminum;  we  shall  use  Zhurkov’s  value  which 
gives,  with  T  =  300  K,f 

a  =  1.4878  x  10 -26/s-  (34) 

We  shall  only  use  eq.  (19a)  from  this  point  on  since  it  is  in  better  agreement  with  the  facts  on  stress 
distribution  than  eq.  (15). 

We  again  use  t  =  1,  c  =  0.25  which  means  we  use  eq.  (31).  Figure  4(a)  shows  the  /?,  s'0  estimates 
plotted  against  time  t,  and  Fig.  4(b)  plots  the  same  estimates  vs  y.  These  graphs  are  similar  in  form 
to  those  shown  in  Fig.  2  for  steel.  The  first  figure  indicates  that  A  is  decreasing  in  a  nonlinear  manner 
with  time;  i.e.  work  hardening  is  a  nonlinear  increasing  function  of  time.  Figure  4(b)  again  shows 
that,  except  in  the  initial  portion,  the  estimates  are  remarkably  close  to  a  straight  line.  Let 
us  examine  this  initial  portion  in  greater  detail. 

The  4a,  =  1  mm  (or  Ay,  =  1/76.2)  in  the  initial  portion  of  Fig.  4(b).  Let  us  use  now 
A  a  ~  0.20  mm.  The  results  are  plotted  in  Fig.  5(a)  for  9  ^  a  13.2  mm.  It  is  clear  that  there  is  a 
departure  from  a  straight  line  in  this  interval.  On  the  other  hand,  if  we  examine  the  interval 
15  ^  a  ^  17.2  mm,  with  A  a  =  0.20  mm,  we  obtain  the  results  shown  in  Fig.  5(b)  where  the  straight 
line  approximation  is  excellent.  We  did  not  encounter  this  departure  from  a  straight  line  in  the 
initial  portion  of  Fig.  2  for  steel.  What  might  cause  this  departure  from  a  straight  line? 

It  is  our  opinion  that  in  the  interval  9  <  a  <  15  mm  the  effect  of  the  load  shedding  program 
is  still  present.  The  reason  we  did  not  see  this  feature  in  Fig.  2  is  that  the  data  considered  started 
at  a0  =  1.0  in.  which  may  have  eliminated  the  effect  of  the  load  shedding  program  that  might  have 
been  present  if  we  used  an  a0  <  1 .0  in.  The  interesting  point  to  observe  is  that  we  now  have  a  method 
for  detecting  when  the  effect  of  load  shedding  ceases. 

We  see  from  Fig.  4(b)  that  if  we  start  at  15  mm.  the  /?,50  estimated  when  plotted  against  y, 
are  essentially  on  a  straight  line.  The  least  squares  estimate  for  this  line  is 

Py  s'o  =  26.20280  (I  -0.72878 y).  (35) 

(c)  IHN  data{  16] 

The  material  is  2024-T3  aluminum  alloy.  The  specimen  is  a  center  crack  panel  400  mm  long, 
70  mm  wide,  and  1  mm  thick.  The  load  is 

F(t)  =  706  +  468  sin  a>t  (lbs)  (36) 

+Actually  T  =  307  K  but  this  is  sufficiently  dose  to  300  K  for  our  purposes. 
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Fig.  5.  estimates  from  VHG  data  plotted  against  yr  \y)  =  0.20/76.2.  Two  /-intervals  are  shown. 

with  co  =  36.8  n  rad/s,  30  replications  were  used.  Data  were  recorded  as  crack  length  to  reach  a 
given  number  of  cycles.  A  load  shedding  program  was  employed  to  produce  the  initial  average 
crack  length  of  aQ  =  7.0302  mm,  /„  =  35  mm.  T  is  assumed  to  be  300  K,  and  we  use  a  from  eq.  (34). 
The  ( t,y )  data  are  available  at  8  points  for  purposes  of  parameter  estimation. 

The  data  analysis[10]  estimated  the  mean  time  to  reach  a  specified  crack  length  which  we  use 
as  in  the  above.  We  use  (19a)  with  e  =  1  and  z  =  0.25.  After  the  first  point,  A_y;  =  2/35. 

Figure  6  plots  the  Pjs'0  estimates  vs  yr  After  the  first  point,  the  estimates  lie  on  a  straight  line 
with  excellent  accuracy.  The  least  square  estimate  of  the  line  is 

fi,Ji  =  26.124396(1  -  0.72293y)  (37) 

where  the  first  point  is  ignored.  The  first  point  does  not  lie  on  the  line  due  to  the  effect  of  the  load 
shedding  program. 

The  material  in  these  data  and  the  Virkler  data  is  2024-T3  aluminum  alloy.  However,  the 
specimen  thicknesses  are  1  mm  and  2.54  mm,  respectively.  This  means  that  the  material  prep¬ 
arations  are  different.  Thus,  the  parameters  in  eq.  (37)  will  be  different  from  those  in  eq.  (35). 

(d)  GD  data 

The  material  is  aluminum  7075-T6  alloy.  The  specimen  is  a  center  crack  panel  320.67  mm  long, 
101.6  mm  wide,  and  3.175  mm  thick;  thus  4  =  50.8  mm.  Three  load  conditions  are  involved: 

Test"  1,  F(t)  =  4099  +  1024  sin  c at  (lb.), 

Test  *2,  F(/)  =  3752  +  1250  sin  tot  (lb.). 

Test*  3,  F(t)  =  2391  +  1024  sin  cot  (lb.),  (38) 


Fig.  6.  0,  estimates  plotted  against  y,  for  the  IHN  data.  A/,  *  2/35  after  first  interval. 
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with  a>  =  20  n  rad/s,  60  replications  were  used  in  each  test.  A  load  shedding  program  was  employed 
to  produce  an  initial  crack  length  u0  =  9  mm.  The  (/.  y)  data  are  available  at  6  points  for  purposes 
of  parameter  estimation. 

These  available  data  analysis  uses  three  Ay/. 

Test'  1,  Ay,  =  0.0360827, 

Test '2,  Ay,  =  0.0393708, 

Test '3,  Ay^  =  0.0328150.  (39) 

The  data  analysis  estimates  the  mean  time  to  reach  6  crack  lengths  starting  at  a0  =  9  mm.  With 
T  =  300  K, 

a  =  1.4878  x  10~26/s 

as  in  the  case  for  data  sets  (b)  and  (c).  Again  we  use  e  =  1,  :  -  0.25,  which  thus  makes  eq.  (31) 
applicable.  Figure  7(a)  shows  the  pjS'0  estimates  plotted  against  y  as  abscissa  for  each  of  the  three 
test  conditions.  Except  for  the  first  point,  the  remainder  can  be  accurately  approximated  by  a 
straight  line  for  each  condition.  A  few  comments  are  in  order. 

First,  the  residual  effect  of  the  load  shedding  program  explains  why  the  first  point  is  below 
the  straight  line  for  each  condition;  the  same  residual  effect  was  noted  for  the  data  sets  (b)  and  (c). 

While  it  is  no  longer  remarkable  to  obtain  these  straight  lines,  it  is  remarkable  that  the  three 
lines  are  almost  parallel.  The  least  square  fit  of  a  straight  line  to  these  estimates,  excluding  the  first 
point,  is 


Test'l: 

-0.71978y), 

)30So  =  25.9862, 

Test'2: 

Pys'0  =  j30ii(l  -  0.72546y), 

P0s'0  =  26.1707, 

Test '3: 

pyn  =  p^{\  —  0.71895y), 

p0s’0  =  25.6854, 

(40) 

the  slopes  being  respectively  -18.7044,  -18.9858,  —18.4665.  These  slopes  are  slightly  different 
but  the  differences  will  not  show  up  on  a  graph.  We  also  note  that  the  coefficient  of  y  in  the 
parentheses  are  very  close  in  value  in  the  three  cases;  the  average  is  0.72140  with  a  maximum 
difference  of  0.56%.  In  view  of  the  fact  that  the  currently  available  data  analysis  only  provides  5 
points  with  which  to  estimate  each  line,  these  numerical  values  cannot  be  viewed  as  final  although 
they  are  interesting. 

Third,  we  know  s’0  k  in  each  case;  thus  we  find 

Test'l:  jg. i  =  9560 psi,  p0  ,=  0.0027 1 8/psi, 

Test'2:  si,  =  9150  psi,  PQ  2  =  0.002860/psi, 

Test '3:  si,  =  6120  psi,  ^0  3  =  0.004197/ psi.  (41) 


Fig.  7.  (a)  p.si,  estimates  plotted  against  y]  for  the  these  load  conditions  of  the  GD  data.  A,v,  variable: 
(b)  /J0  estimates  plotted  against  si,  for  these  load  conditions. 
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The  /}„  k  estimates  are  plotted  against  s, in  Fig.  7(b).  We  note  these  estimates  lie  close  to  a  straight 
line,  although  the  points  are  not  well  spaced.  The  least  square  linear  fit  obtained  from  eq.  (41)  is 

P0  =  0.006838  -  4.32535  x  10''s0.  (42) 

Thus.  (42)  suggest  that  given  any  S'n  we  can  calculate  f}0  for  this  s'n.  Since  all  coefficients  of  y  in 
(39)  are  close  to  0.72140,  this  further  suggests  that  prediction  of  mean  crack  growth  at  s,j  different 
from  those  tested  might  be  possible.  More  will  be  said  about  this  below. 

Finally,  we  observe  in  eq.  (41)  that  the  /?„  are  different  for  different  s'0.  Since 


(43) 


the  results  in  eq.  (41)  say  that  A0  depends  on  the  loading  through  T,',  This  means  that  the  work 
hardening  process  depends  on  the  loading  as  well  as  on  the  dimensionless  crack  length  y  as 
indicated  by  eq.  (40).  We  thus  see  that  /?,  5 <j  depends  on  the  two  fixed  material  parameters  /?„  or 
A0  and  the  coefficient  of  y  in  eq.  (40),  of  which  /?0  depends  strongly  on  the  loading  through 
whereas  the  coefficient  of  y  does  not  seem  to  do  so. 

The  GD  data  are  the  most  interesting  of  the  four  sets  (a),  (b),  (c)  and  (d)  considered  in  this 
section,  because  of  the  three  sets  of  loads  employed.  However,  a  data  analysis  with  more  than  15 
points  in  each  of  the  three  crack  length  intervals  is  needed  to  improve  the  accuracy  of  the  values 
given  in  eqs  (40)  and  (41). 

Let  us  next  discuss  what  can  be  said  when  the  four  data  sets  are  considered  together. 


6.  DISCUSSION 

The  consistency  of  the  character  of  the  estimates  of  /?,  from  four  data  sets  involving  different 
metals,  load  means  and  amplitudes,  and  laboratories  is  certainly  remarkable.  We  have  already 
discussed  in  the  previous  section  some  of  the  physical  inferences  that  can  be  drawn  from  these 
consistent  results  on  mean  behavior  in  metals  where  geometry  of  the  crack  development  is 
essentially  the  same.  There  are  many  points  that  merit  further  discussion.  However,  in  this  initial 
paper,  we  only  can  take  up  a  few. 

Equations  (8)  and  (9)  are  the  starting  point  for  this  study  at  the  macro  or  engineering  level. 
Equation  (9)  is  a  force  balance  equation  and  requires  no  further  comment  here.  Equation  (8)  says 
that  two  ways  of  evaluating  the  macro  mean  rate  of  bond  breaking  are  equal,  and  relates  a(t)  to  the 
bond  breaking  process  in  the  unbroken  portion  of  the  cross-section.  It  is  important  to  realize  that  the 
left-hand  side  of  eq.  (8)  implies  that  all  the  bond  breaking  takes  place  at  the  crack  tip  or  at  least 
in  a  small  region  just  ahead  of  the  crack  tip.  Obviously,  the  right  hand  side  must  say  the  same  thing. 

We  approximated /in  eq.  (3)  by  means  of  the  normal  stress  s(t,  ?)•  Stress  is  an  average  concept. 
Love[16]  has  pointed  out  that  it  is  a  “macroscopic  or  large-scale"  concept,  which  means  we  are 
using  a  smooth  version  of /.  While  this  is  the  simplest  method  of  approximating/ for  the  geometry 
and  loading  under  consideration,  it  may  not  be  the  only  possibility  that  merits  attention[17]. 
Further,  with  different  specimens  geometry  and  method  of  applying  the  load,  the  normal  stress 
s(t,  Z)  will  not  be  appropriate,  and  we  leave  this  point  to  another  time.  Let  us  next  focus  attention 
on  the  normal  stress. 

The  choice  of  s(t,  ,;)  in  eqs  (8)  and  (9)  is  basic  to  what  follows.  As  noted  earlier.  s(t,  ;)  always 
occurs  under  an  integral  sign  with  respect  to  and  this  means  that  the  broad  features  of  its  form 
as  a  function  of  J;  are  important  but  not  the  details.  When  a  macro  crack  is  present,  a  uniform 
normal  stress  is  not  realistic.  Equation  (17)  does  contain  the  broad  features  that  must 

possess,  and  this  is  the  simplest  choice  that  can  be  made  possessing  these  features  (see  ref  [18]  for 
a  different  approach).  However,  this  choice  requires  the  addition  of  the  two  new  quantities  A  or 
c)  and  r7  when  eq.  (17a)  is  employed.  At  this  time,  there  are  no  experimental  data  to  suggest  how 
these  two  quantities  behave.  Qualitatively,  we  know  that  the  plastic  zone  increases  in  length  and 
work  hardening  increases  as  time  or  crack  length  increases  for  the  materials  considered.  However, 
we  have  no  quantitative  information  on  how  v,  <5.  and  r„  jointly  evolve  with  time  or  if  they  may 
be  regarded  as  a  state  vector(lO).  To  circumvent  this  current  lack  of  information  and  get  on  with 
a  preliminary  parameter  estimation  study,  we  assume  the  constraint  eq.  (29),  which  is  at  least 
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physically  plausible.  Clearly,  it  is  important  to  obtain  experimental  data  on  how  y.  S  and  r„  jointly 
evolve  with  time  in  order  to  refine  and.  or  expand  the  scope  of  this  model. 

We  employ  the  same  s(t.  -)  for  center  crack  panel  and  WOL  specimens  under  tensile  loading. 
This  choice  can  be  questioned.  However,  we  see  from  the  data  analysis  with  £  =  1  and  r  =  0.25. 
that  the  high  stress  region  produces  the  overwhelming  fraction  of  bond  breaking  that  takes  place. 
This  would  still  be  the  case  if  £  ^  0.20.  Thus,  the  precise  form  of  j(t,  in  the  low  stress  region, 
which  will  differ  between  the  specimen  types,  is  not  important. 

The  crack  growth  process  is  more  complex  than  failure  under  constant  stress[3].  To  see  this, 
let  us  write 


/U;  =  j?05o(i  -Tv).  (44) 

which  conforms  to  all  cases  considered.  The  substitution  of  eq.  (44)  into  the  first  term  of  eq.  (30) 
gives 

Unroll  -<)i I  I 

Ay,  =  -j  y,  At, e  1  .  (45) 

which  implies  e  >  0.20.  Instead  of  the  two  constant  material  parameters  a  and  /(  employed  in  refs 
[2,  3],  we  now  have  three,  namely,  a.  ji0  and  c.  The  parameter  a  depends  on  the  material  but  not 
on  the  load;  /?0  depends  on  the  material  and  the  load;  and  c  depends  on  the  material  but  not  strongly 
on  the  load.  Thus,  the  parameter  estimation  problem  and  the  model  defined  by  eq.  (45),  which  uses 
eq.  (29),  is  more  complex  than  the  models  considered  by  previous  investigators.  This  fact  has  been 
established  with  the  aid  of  four  high  replication  sets  of  FCG  data. 

The  presence  of  T  in  an  explicit  manner,  in  r  and  /?0,  is  important.  Its  presence  indicates  that 
it  should  be  measured  during  tests.  Further,  since  relatively  small  changes  in  T  can  produce 
substantial  changes  in  fatigue  life,  lack  of  knowledge  of  how  T  may  be  changing  can  introduce  a 
source  of  life  variability  that  may  produce  incomprehensible  results  as  well  as  unwanted 
consequences.  Thus,  having  T  explicitly  present,  is  a  source  of  strength  for  the  current  model. 

The  purpose  of  the  load  shedding  program  in  FCG  tests  is  to  eliminate  the  transient  material 
behavior  that  is  induced  when  the  initial  crack  length  is  rapidly  produced.  Consider  the  upper  graph 
in  Fig.  5  which  refers  to  the  15  to  17  mm  crack  length  interval.  The  fylo  estimates  essentially  lie 
on  a  straight  line,  indicating  in  our  view,  that  the  rate  of  work  hardening  is  constant.  We  interpret 
this  to  mean  material  behavior  has  reached  a  steady  state.  The  lower  graph,  which  is  for  the  9  to 
13  mm  interval,  shows  an  initial  departure  from  steady  state  material  behavior.  We  believe  this 
effect  is  due  to  the  lingering  effect  of  the  load  shedding  program.  This  appears  to  be  the  first  model 
that  can.  from  a  data  analysis,  detect  this  lingering  effect.  Hence,  this  is  a  remarkable  feature  of 
the  model.  There  is  another  point  worthy  of  comment. 

Consider  the  lower  graph  in  Fig.  4,  which  is  for  the  9  to  49  mm  crack  length  interval.  The 
straight  line  character  of  the  estimates  persists  to  the  largest  a -values  employed,  where  the  crack 
tip  is  coming  close  to  the  panel  edge.  The  same  feature  is  observed  in  Figs  2,  6  and  7.  It  is  reasonable 
to  expect  that  the  close  proximity  of  the  crack  tip  to  the  specimen  edge  would  disturb  the  constant 
rate  of  work  hardening  in  view  of  the  fact  that  the  lingering  effect  of  load  shedding  does  change 
the  rate.  Thus,  it  is  a  surprise  that  there  is  no  change  in  the  rate  due  to  the  edge  effect,  even  in 
Fig.  2  where  the  crack  tip  comes  to  within  0.55  in.  of  the  specimen  edge. 

Numerous  approximations  were  made  in  obtaining  eqs  (27)  and  (31).  For  example,  the  varying 
y-value  on  a  Ay,  interval  was  replaced  by  the  constant  y y,  sin  r  was  replaced  by  a  square  wave 
function;  and  because  of  the  numerical  values  of  the  parameter  being  estimated  and  our  choice 
£  =  1,  certain  terms  became  negligible  and  were  dropped.  In  a  preliminary  study,  these  approxi¬ 
mations  are  acceptable.  However,  now  that  we  have  found  that  the  model  is  interesting,  it  becomes 
necessary  to  refine  our  parameter  estimation  study  by  eliminating  these  approximations,  and  this 
will  be  done  in  a  future  study. 

The  model  presented  is  for  mean  behavior.  However,  we  mentioned  earlier  that  if  variability 
is  to  be  included,  we  must  re-examine  our  evaluation  of  the  term  A /,  although  this  term  may  not 
be  the  only  source  of  variability.  What  makes  this  term  of  interest  is  that  /  and  A  vary  from  bond 
to  bond  due  to  variability  in  the  atomic  structure  of  the  material.  Further,  due  to  the  polycrystalline 
nature  of  the  material,  different  crystals  and  grain  boundaries  will  be  stronger  than  others  providing 
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possible  correlation  in  /.-values  over  intervals  of  varying  lengths.  These  possibilities  are  some  that 
will  be  of  interest  when  variability  is  taken  into  account,  which  we  also  leave  to  a  future  paper. 


7.  CONCLUSION 

The  model  presented  for  mean  fatigue  crack  growth  is  based  on  fundamental  results  provided 
by  physical  theory.  The  starting  point  for  the  derivation  of  the  model  is  the  basic  expression  for 
the  mean  rate  of  single  atomic  bond  breaking  provided  by  the  kinetic  theory  of  thermally  activated 
bond  fracture.  Based  on  this  expression,  which  applies  at  the  micro  level,  a  macro  model  is  derived 
for  mean  fatigue  crack  growth,  as  a  function  of  time,  in  a  center  crack  panel  and  WOL  test 
specimen  under  (simple  harmonically  varying)  tensile  load.  The  parameters  in  the  model  have 
physical  meaning  due  to  either  the  time  varying  load  or  material  properties.  In  this  paper,  only 
the  simplest  version  is  considered. 

The  ideas  employed  here  differ  in  a  fundamental  manner  from  those  employed  in  fracture 
mechanics.  In  this  paper,  basic  ideas  from  physics  are  brought  in  at  the  atomic  level  at  the  outset. 
This  is  different  from  the  traditional  approach  and  offers  a  basic  approach  that  has  not  been 
previously  studied.  The  model  is  not  phenomenological. 

The  most  significant  point  of  the  macro  equations  and  parameter  estimates  that  have  been 
developed  is  their  engineering  importance.  The  most  important  feature  for  the  engineer  is  that  for 
a  specific  component  only  two  or  at  most  three  tests  are  required  to  predict  mean  fatigue  growth 
for  an  arbitrarily  varying  simple  harmonic  tensile  load. 

We  see  from  eq.  (45)  that  there  are  three  parameters,  namely,  a,  /?0  and  c  that  must  be  estimated 
in  order  to  use  this  equation  for  prediction  of  mean  fatigue  life  behavior.  The  parameter  a  can  be 
estimated  from  tabulated  data.  The  parameters  /?0  and  c  must  be  estimated  from  tests.  While  P0 
depends  on  ja  and  As  through  s$,  c  does  not  seem  to  depend  on  these  quantities.  We  observe  from 
the  GD-data  analysis  that  two  or  at  most  three  suitably  chosen  tests  are  required  to  estimate  the 
dependence  of  /?„  on  the  load;  c  is  then  estimated  as  the  average  of  the  test  estimates.  Once  this 
is  done,  we  can,  on  selecting  suitable  values  for  t  and  r,  employ  eq.  (45)  for  predictive  purposes 
for  an  arbitrarily  varying  simple  harmonic  tensile  load  at  any  temperature.  This  predictive 
capability  based  on  limited  testing  is  a  significant  feature  of  engineering  importance. 

Many  interesting  points  concerning  the  mean  model  merit  exploration  and  discussion.  Further, 
the  variability  inherent  in  the  phenomenon  must  be  brought  into  a  model.  We  leave  these  matters 
for  future  papers. 
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APPENDIX 

In  Fife.  Al,  we  plot  points  from  EN  vs  a  and  from  EA  vs  n.  These  points  come  from  the  data  analysis  given  in  [10] 
We  observe  that  for  all  practical  purposes  the  curves  these  sets  of  points  lie  on  are  the  same  curve.  This  justifies  the  use 
of  EN  vs  a  in  the  above  data  analysis  instead  of  the  EA  vs  n  curve. 


Fig.  AI.  A  comparison  of  EN  vs  a  points  and  EA  vs  n  points  for  VHG  data. 
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Abstract 


A  study  is  continued  of  a  macro  model  of  FCG  that  is  based  upon  a  micro  result  from 
reaction  rate  theory.  The  previously  presented  deterministic  (mean)  model  is  examined  further, 
and  a  probabilistic  model  is  introduced.  Analysis  of  fatigue  crack  growth  (FCG)  data  for 
periodic  tensile  loads  reveals  that  three  material  parameters  are  required.  The  model  also 
explicitly  contains  the  load  and  the  temperature  T.  It  is  shown  that  T  and  the  form  of  the 
periodic  load  are  important  quantities  of  which  the  former  has  interesting  implications  for 
accelerated  testing.  The  minimum  amount  of  testing  required  for  predictive  purposes  is  also 
discussed. 
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1.  Introduction 

We  present  in  [1]'  the  derivation  of  a  deterministic  macro  model  of  fatigue  crack  growth 
(FCG)  based  on  a  micro  result  from  reaction  rate  theory  along  with  an  initial  study  of  some  of  its 
properties.  This  presentation  provides  a  brief  survey  of  some  of  the  interesting  features  of  the 
FCG  process  and  its  modeling  by  this  approach. 

The  purpose  of  the  present  paper  is  threefold:  First,  further  studies  of  the  deterministic 
(mean)  model  for  the  evolution  of  the  dimensionless  crack  length  are  presented  and  discussed. 
We  use  for  this  purpose  a  data  analysis  of  six  sets  of  FCG  in  which  the  specimens  are  under 
simple  harmonically  varying  tensile  load.  The  studies  indicate  the  importance  of  the  model 
having  three  material  parameters,  one  of  which  depends  on  the  load.  The  model’s  explicit 
dependence  on  the  temperature  T  and  the  wave  form  of  the  periodic  load  also  indicates  the 
importance  of  these  two  quantities  in  the  evolution  of  the  crack  length. 

Second,  and  initial  study  is  made  of  how  the  ever  present  variability  in  FCG  might  be 
incorporated  into  the  macro  mean  or  deterministic  model  by  means  of  random  material  behavior. 
A  data  analysis  of  the  variability  present  in  the  data  sets  is  employed  to  estimate  the  distributions 
of  random  material  behavior,  such  distributions  then  producing  the  distributions  of  the  time  to 
reach  a  given  crack  length.  The  closeness  of  the  standardized  distributions  of  the  random 
material  behavior  is  discussed  and  suggests  that  there  might  be  an  underlying  micro  random 
material  behavior  that  merits  study. 

Finally,  we  consider  the  minimum  amount  of  testing  required  to  estimate  the  macro  model 
parameters  for  predictive  purposes,  and  the  most  expeditious  method  for  accelerating  tests. 
Since  only  one  material  parameter  depends  on  the  load,  it  turns  out  that  two  or  three  tests  at 
suitably  chosen  loads  are  needed  for  parameter  estimation;  then  it  is  possible  to  predict  the  time 

*  Hereinafter  referred  to  as  Part  I  or  simply  as  I. 
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to  reach  a  given  crack  length  over  a  wide  range  of  loads.  A  standard  method  for  accelerating 
fatigue  and  FCG  tests  is  to  increase  the  loads,  which  introduces  the  question  of  the  effect  of 
changed  load  interaction.  The  model  indicates  that  testing  time  can  be  shortened  by  the  simple 
expedient  of  increasing  the  temperature  T  with  the  same  loads,  which  eliminates  the  load 
interaction  question. 

We  begin  by  giving  a  short  review  of  the  derivation  of  the  deterministic  model. 

We  close  with  a  conclusion  section  that  summarizes  the  results  of  the  study. 
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2.  Brief  Review  of  Deterministic  Model 

The  basic  specimen  geometry  under  consideration  is  a  thin  center  crack  panel  under  axial 
tensile  load  F(t)  (See  Figure  1  of  I).  We  start  from  the  macro  equations  (81)  and  (91): 

—  =[’°aeps<^d^.  (1) 

dt  Ja 

F(t)  =  2wfS(t,')d^,  (2) 

where 


a  =  1/2  crack  length, 
l0  =  1/2  panel  width, 
w  =  panel  thickness, 
q  =  distance  from  panel  center  line  to  a 

generic  point  in  unbroken  portion  of  crack  plane, 

AG* 

a  =  — -e  kT  ,  B  = - ,  T  =  absolute  temperature  , 

h  vkT 

k  =  Boltzmann's  constant ,  h  =  Planck's  constant , 

a.  =  average  distance  over  which  force  acts  in  bond  breaking  , 


(3) 


v  =  mean  number  of  atomic  bonds  per  unit  area  . 

s(t,  q)  is  the  normal  stress  on  the  unbroken  portion  of  the  crack  plane  at  time  t  and  location  it 
is  a  bounded  function.  AG4-  is  a  material  parameter  that  can  usually  be  estimated  from  existing 
information.  The  growth  control  parameter  (3  or  X/\)  depends  on  the  material  and  how  it  is 
formed  into  components;  it  must  be  estimated  from  data. 
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Neglecting  dynamical  effects  in  the  specimen,  we  can  assume  (see  ( 1 61)) 

s(t,  q)  =  s(c)g(t) , 


(4) 


where  the  stress  factor  s(*)  only  depends  on  q  and  g(t)  is  a  dimensionless  multiplicative  time 
factor  depending  on  t  and  crack  length. 

The  normal  stress  is  high  in  the  region  of  irreversible  bond  breaking  just  ahead  of  the 
crack  tip  and  then  decreases  to  an  elastic  stress  as  we  move  away  from  the  crack  tip  in  the  crack 
plane.  A  realistic  assumption  for  s(q)  is 


s(q)  =  cp  ,  a  <  q  <  a  +  A  , 

(5) 

=  ae(<  ap),  a  +  A  <^</0  , 

(see  (171)).  Other  assumptions  are  possible  (see  (17bl)),  but  (5)  is  the  simplest  in  view  of  the 
fact  that  (1)  and  (2)  are  not  overly  sensitive  to  the  precise  form  of  s(q). 

The  substitution  of  (4)  and  (5)  into  (1)  and  (2)  yields 


i(t)  = 


s0(t)/ae 


1  -  y  -  5  +  ra5 

Pr3s0(t) 


dv 

dt 


=  a 


Se 


l-v-0  +  r-5 


+  (1  -  y  -  8)e 


gso(0 
l-v-5  +  ra5  ! 


1 


(6) 

(7) 


where  s0(t)  =  F(t)/2w/0  is  the  far  field  axial  tensile  stress,  y  =  a JlQ  is  the  dimensionless  crack 
length,  5  =  A//0  is  the  dimensionless  length  of  the  high  stress  region  in  the  crack  plane,  and 
ap/ce  =  ra  =  1  -r  €  is  the  ratio  of  the  high  stress  to  the  elastic  stress  and  is  greater  than  one.  Of 
course,  1  -  y  -  5  +  raS  >  0  and  0  <  y  <  1. 

Given  a,  (3  and  sQ(t);  (7)  must  be  integrated  to  yield  y  as  a  function  of  time.  However,  for 
the  specific  form  of  s0(t)  we  shall  assume,  a  difference  form  of  (7)  is  satisfactory  for  our 


purposes. 
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Let  the  far  field  axial  stress  s0(t)  acting  on  the  specimen  be 

s0  (t)  =  ?0  +  As  sin  cot ,  or  (8) 

=  ?0  +  Ash(t) , 

where  r  =  cot,  h(r)  is  continuous  and  periodic  with  period  2;t,  |  h(x)  |  <  1,  it  is  odd  with  respect  to 
x,  so  its  mean  value  is  zero,  and  ?0  and  As  >  0.  Let  us  discretize  the  crack  length 
a0  <  ai  <  •  •  •  <  a n  <  /0,  and  time  0  <  q  <  •  •  •  <  tn,  where  tj  is  the  mean  time  to  reach  aj. 
Further,  let  us  set 


Vj=aj//0  . 

yJ  =  (yJ-i+yj)/2  . 

Ayj  =  y j  -  Vj_i  , 


(9) 


~  (j-i  * 

where  At,  in  the  time  required  for  y  to  increase  by  Ay^.  The  difference  form  of  (7)  that  we 
employ  here  is 


AVj  =  ocAtj 


1  f2;t  Vjfq^shCT) 

2k  •'o 


dx  +  (1  -yj  -5j)ev,s° 


1  r2^  VjAsh(t) 

2k  o  C 


dt 


(10) 


where  =  Pj/(1  -  yj  -  8j  +  ra5j),  co  converts  cycles  An;  to  time  Atj. 

The  parameters  ra  and  5  describe  how  the  assumed  form  of  the  normal  stress  on  the  crack 
plane  evolves  with  time.  It  is  possible  that  y,  r0,  and  5  form  a  state  vector  if  we  have  differential 
equations  describing  the  joint  evolution  of  y,  ra,  and  5.  This  would  imply  that  y  is  history 
dependent.  However,  at  this  time,  we  have  no  such  evolutionary  equations  for  these  quantities. 
Therefore,  to  proceed  with  a  study  of  (10)  in  conjunction  with  data  we  introduce  the  constraints 
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ra  =  1  -f  e  ,  e  =  constant  >  0  , 


5j  =  zv,  ,  z  =  constant  with  0  <  z  <  1  . 


(11) 


The  first  constraint  says  that  as  ae  increases  with  increasing  y  ap  also  increases  at  the  same 
rate.  The  second  constraint  says  that  the  length  of  the  high  stress  region  of  irreversible  bond 
breaking  in  the  crack  plane  increases  as  y  increases.  These  statements  are  consistent  with  the 
qualitative  behavior  of  a  work  hardening  material.  Eq.  (11)  are  among  the  simplest  that  can  be 
introduced  but  can  be  changed  if  necessary.  More  will  be  said  about  (1 1)  when  confronting  data. 

If  e  >  .20,  parameter  estimation  studies  reported  in  I  and  reported  here  indicate  that  the 
second  term  in  the  square  bracket  of  (10)  is  negligibly  small  in  comparison  with  the  first  term. 
Thus,  (10),  when  combined  with  '  1  *)  can  be  written  as 


2k  vj  =  azVjAtje 


vj(1-»-e)s0 


J_f2Vi(1+e)Ash(t)dX 

2jc  ■'o 


(12) 


We  coitoider  two  forms  of  h(x),  namely 


h(t)  =  sin  x , 


h(x)  =  —  ,  0  <  x  <  rt/2  , 

K 


2x 


-  n  -  —  ,  k/2  <  r  <  3n/2  , 

K 


(13) 


=  -  4  +  —  ?  37T/2  <  x  <  2k  , 

K 

The  second  form  is  a  "saw-tooth'’  version  of  the  first  each  having  amplitude  1.  First,  consider 
the  first  of  '  IT);  then 


where  I0(x)  is  the  modified  Bessel  function  of  order  zero.  For  the  data  studied,  x  is  large;  thus 
w'e  can  employ  the  approximation  [2] 


ex  (1  +£)PjAs 

IoW-  -  ,  X= - — 

M2tzx  1  —  (1  —  z  £)vj 

Hence,  we  can  approximate  (12)  by  means  of 

Avj  =  azyjAtj(2^UjAs,)-1/ieUi  ,  or  Atj  =  (azyj)-1  Avj(2:cujAs  )I/ie~Uj 

when  As  is  not  too  small  (e.g.  x  >  10),  and 

(1  +  e)p,s0 ' 

Uj  = - —  , 

1  -  (1  -  z  e)vj 

_  /  __ 

Sq  —  Sq  +  As  , 

As  =  As/T0  . 


(15) 


(16) 


(17) 


Next,  consider  the  second  form  of  h(t)  given  in  (13).  Then,  under  the  same  type  of 
approximation  we  find 

azv.Atj  ,  .  u 

Avj  =  — —  (UjAs  )-1e  !  ,  (18) 

where  (17)  still  holds.  We  assume  in  Pan  I  that  h(x)  =  rectangular  wave  form  with  the  area 
under  the  half-w-ave  equal  to  the  area  under  sin  t  for  0  <  r  <  k.  The  equation  corresponding  to 
(16)  and  (18)  then  becomes 


azyjAtj  u 
Ayj  = - — ' 


(19) 


which  is  (311)  if  z  =  1/4,  and  s0  =  s0  +  — As,  e  =  1.  In  (18)  and  (19),  As  is  not  too  small 

K 
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Equations  (16),  (IS),  and  (19)  are  based  on  the  step  distribution  for  the  stress  s(q)  given  by 
(5).  Another  possible  choice  is  given  by  (17bl);  still  other  reasonable  possibilities  exist. 
However,  at  this  point  in  the  study,  use  of  (5)  is  appropriate  because  of  its  simplicity. 

The  material  parameters  in  (16),  (18),  and  (19)  are  a  and  [3;  e  and  z  are  geometrical 
parameters  connected  with  (5).  We  have  already  stated  that  the  parameter  a  (or  AG+)  can  be 
estimated  from  known  results  (See  Zhurkov  [3]).  However,  the  pj  must  be  estimated  from  data. 
The  selection  of  values  for  €  and  z  will  be  discussed  below. 

Let  us  consider  how  pj70  will  be  estimated  from  data,  using  (16).  From  the  data,  we  can 
find  for  each  Avj  Atj  =  E(ATj },  where  ATj  is  a  random  variable  representing  the  time  for  the 
crack  to  increase  by  Avj  in  dimensionless  length.  Then,  write  (16)  as 


AVjaST  _  eUi 
azyjAtj  a/ujAs 


(20) 


On  taking  logarithms,  we  obtain 


1  Av;  i 

—log  2k  +  iOg — - - =  uj  -  —log  u_jAs  . 

azyjAtj 


(21) 


Set 


uj  =  bj  +  dj 

It  follows  from  (21)  that 

dj  =  ~-log  bjAs'  +  -i-log  (1  +  ~) 

_  Dj 

and  since  dj/bj  is  small  for  the  data  considered 

dj  -  (1-^-r1  {log  bjAs' 


(22) 


(23) 


(24) 
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Thus,  we  have  an  estimate  for  u, .  We  then  find  an  estimate  for  P,s0  by  employing  the  first  of 
(17) 


PjS0 


1  -  (1  -  Z  £)Vj 
1+8 


(25) 


once  z  and  €  are  selected.  Similar  results  can  be  obtained  for  (18).  We  easily  find  from  (19) 
that 


2Av; 

Uj  =  log — -li-  ,  (26) 

azvjAtj 

which  gives  on  estimate  for  Uj  in  this  case;  estimates  for  (3j?0  again  follow  from  (25). 

Eq.  (7)  is  the  basic  contribution  in  this  paper  (and  I).  It  contains  the  significant  quantities 
of  engineering  interest.  Comparison  with  two  other  models  of  FCG  based  on  reaction  rate 
theory  will  be  commented  upon  in  Section  4. 
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3.  Parameter  Estimates  for  Deterministic  Model 

We  employ  tl.e  same  experimental  data  on  FCG  that  is  considered  in  I,  namely 

a)  A  STM  E24.04.0 1  Round  Robin  Data  [4] 

The  material  is  lONiSColMo  steel,  24  replications,  co  =  1071  rad/sec. 

b)  VHG  Data  [5] 

The  material  is  2024-T3  al  alloy,  68  replications,  co  =  407:  rad/sec. 

c)  EHN  Data  [6] 

The  material  is  2024-T3  al  alloy,  30  replications,  co  =  36. St:  rad/sec. 

d)  GD  Data  [7],  3  load  conditions,  co  =  207:  rad/sec. 

The  material  is  7075-T6  al  alloy,  60  replications  for  each  load  condition. 

Let  tjk  be  the  time  to  reach  crack  length  aj  for  sample  k  in  one  experiment,  where  k  =  1,  •  •  • ,  n. 
Then,  with 

Atj,k  =  tj.k  ~  >  (27) 

we  employ 

Atj  =  |iAt,k  (28) 

n  k=l 

as  the  (mean)  time  required  for  the  crack  length  to  increase  from  arl  to  3j  or  to  increase  by 
amount  AUj  =  aj  -  aj_i .  Thus,  for  each  load  condition,  we  have  Atj  and  Avj  at  a  number  of  crack 
length. 

We  employ  (16)  for  parameter  estimation  since  in  all  experiments  h(x)  =  sin  x.  In  b),  the 
controlled  ambient  air  temperature  is  297°  K  and  we  assume  a  10°  C  rise  in  specimen  temperature 
due  to  material  hysteresis  at  20  Hz,  producing  an  assumed  steady  state  T  =  307° K.  In  a),  c),  and 
d),  we  assume  as  a  reasonable  value  the  temperature  T  =  300°K. 
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The  AG*  value  in  each  experiment  is  given  in  Pan  I.  Thus,  for  each  experiment,  we  have 
an  estimate  for  a  from  the  formula  given  in  (3),  and  needed  in  (16). 

The  quantities  ?0  and  As  are  obtained  from  the  load  condition  and  specimen  geometry 
for  each  experiment  by  means  of  the  second  and  third  of  (17);  these  quantities  are  also  needed  in 
(16). 

/ 

It  follows  that  given  the  (Atj,  AVj)  (16)  provides  the  estimates  for  the  parameters  pj?0,  for 
the  selected  values  of  z  and  e;  then,  we  get  the  pj  from  these  estimates. 

WOL  specimens  are  employed  in  data  set  a),  whereas  in  the  other  data  sets  the  specimens 
are  center  crack  panels.  Since  significant  bond  breaking  only  occurs  in  the  high  stress  region,  we 
will  employ  the  results  in  Section  2  for  this  case  even  though  rc  and  ae  must  be  viewed  in  a 
different  manner  than  for  center  crack  panel  specimens. 

The  general  nature  of  the  parameter  estimates  Pj?0  found  here  is  of  the  same  as  found  in 
Part  I,  where  we  assumed  h(x)  =  rectangular  wave  form.  That  is,  in  the  data  sets  of  b),  c),  and  d), 
there  is  an  initial  deviation  of  the  (3j?0  from  a  straight  line  with  negative  slope,  which  we 
ascribed  in  I  to  transient  material  behavior,  no  such  deviation  from  a  straight  line  occurs  in  data 
set  a).  Thus,  only  a  few  of  the  graphs  presented  in  I  will  be  repeated  here.  Instead,  we  confine 
attention  to  the  expression  for  the  straight  line  approximation  to  the  estimates. 
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a)  E24.04.01  Data 

Figure  1  shows  the  parameter  estimates  (3j70  from  Part  I  where  h(x)  =  rectangular  and  the 
estimates  from  (18),  both  with  e  =  1  and  z  =  1/4.  The  crack  length  increases  from  1  in.  to  2  in. 
and,  with  Aaj  =  .05  in.,  this  produces  twenty  aj;  /0  =  2.55  in.  and  we  assume  T  =  300°K.  The 
straight  line  approximations  differ  by  a  significant  amount.  The  mean  square  (ms)  approximate 
lines  obtained  with  the  use  of  (18)  for  estimating  the  [3j?0  is 

pys0  '  =  59.534265  (1  -  .723460  y) ,  (29) 

where  ?0  =  1725.49  psi  and  As  =  .450000.  The  corresponding  equation  from  I  is 

Py?0 '  =  59.2759  (1  -  .73625  y)  (30) 

As  noted  above,  there  is  no  deviation  of  the  estimates  from  a  straight  line  when  starting  at 
aj,  =  1  in.  and  ending  at  an  =  2  in.  even  though  at  the  end  point  the  crack  tip  is  close  to  the 
specimen  boundary. 
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b)  VHG  Data 

The  straight  portion  of  the  estimates  is  from  crack  length  16  mm.  -  37  mm.  (Aa;  =  1  mm.), 
37  mm.  -  44.2  mm.  (Aaj  =  1.2  mm.),  44.2  mm  .  -  49  mm.  (Aaj  =  1.6  mm.),  and  49mm-49.8mm 
(Aaj  =  .8  mm.),  giving  31  values  for  aJ5 /0  =  76.2  mm.  so  that  yj  =  aj/76.2  and 
AVj  =  Aaj/76.2,  T  =  307°K.  The  ms.  approximate  line  is 

pyT0  '  =  25.952791  (1  -  .725106  y) ,  (31) 

where  ?0  =  8,766.67  psi.  As  =  .401 141,  e  =  1,  and  z  =  1/4.  The  initial  portion  from  9  mm.  -  16 
mm.  provides  estimates  that  deviate  from  a  straight  line  in  the  same  manner  as  shown  in 
Figure  4  and  5  in  Part  I.  The  equation  from  I  corresponding  to  (31),  when  h(x)  has  a  rectangular 
wave  form,  is 


Pys0  '  =  26.20280  (1  -  ,72878y)  (32) 

Thus,  there  is  no  reason  for  graphing  the  estimates. 

In  the  first  crack  length  interval  16  mm.  -  37  mm.,  we  use  Aaj  =  1mm..  However,  the  data 
in  that  interval  are  given  at  Aaj  =  .2  mm..  This  provides  an  opportunity  for  checking  on  whether 
or  not  use  of  smaller  Ayj  produces  a  significant  difference  in  the  ms.  approximate  line.  On 
comparing  the  ms.  approximate  lines  in  the  interval  16  mm.  -  36  mm.,  using  Aaj  =  1  mm.  and  .2 
mm.,  we  find  that  the  numerical  values  in  the  ms.  approximate  lines  differ  by  a  maximum  of 
.04%.  It  thus  appears  that  using  smaller  Avj  does  not  produce  a  significant  difference  in  the 
parameter  estimates.  However,  a  bigger  difference  in  parameter  values  appears  when  the  Aaj 
were  changed  from  1.0  mm.  to  2  mm.  Thus,  the  magnitude  of  Aaj  should  not  be  too  large. 
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c)  IHN  Data 

The  straight  portion  of  the  estimated  pj?0  is  for  aj  =  9  mm.,  11  mm.,  •  *  • ,  21  mm.,  with 
Aaj  =  2  mm.,  giving  7  points.  /0  =  35  mm.  and  T  =  300°K.  The  ms.  approximate  straight  line  is 

s0  '  =  26.667678  ( 1  -  .723869  y) ,  (33) 

where  I0  =  10,814  psi.  As  =  .397633,  6  =  1,  z  =  1/4.  The  estimates  are  obtained  using  (16)  and 
the  (Ayj,  At;)  are  obtained  from  the  data.  The  estimates  at  aj  =7.5  mm.  deviates  from  the 
straight  line  (33)  as  is  the  case  in  I  (See  Figure  6  of  I). 
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d)  GD  Data 

We  assume  e  =  1  and  z  =  1/4.  In  Part  I  (see  Figure  7  I),  we  employ  six  aj.  Here  we  use 
24  a,  covering  the  same  range  of  crack  lengths  for  the  three  load  conditions  as  in 

I.  /0  =  50.8  mm.  and  T  =  300°K.  The  24  (3j?0  estimates  are  obtained  using  (i6).  Figure  2 

shows  the  [3j?0  estimates  for  the  three  load  conditions. 

The  first  five  estimates  show  transient  material  behavior  in  each  case.  The  ms. 
approximate  straight  lines  for  the  straight  position  are,  using  19  aj  in  each  case, 

GD*  1  :  PySo '  =  26.712158  (1  -  .722287  y) , 

GD*2  :  (3y?0  '  =  26.954427  (1  -  .727481  y)  ,  (34) 

GD*3  :  (3yS0  '  =  26.518554  (1  -  .722104  y) , 

where 

GD#1  :?0  =  10,246  psi  ,  As  =  .199883, 

GD#2  :s0' =  10,004  psi  ,  As' =  .249900  ,  (35) 

GD*3  :  s0'  =  6,830  psi  ,  As' =  .299854  . 

Equations  (34)  and  (35)  produce 

GD#  1  :  p0  =  .00260708  , 

GD*2  :  p0  =  .00269436  ,  (36) 

GD*3  :  p0  =  .00388265  . 

The  form  of  (34)  is 

Pys0'  =  p0s0'(l-cy)  (37) 

The  average  value  of  c  is  .723957,  with  a  maximum  deviation  of  .4867%.  It  appears  that  c  does 
not  depend  strongly  on  the  load  condition.  This  same  observation  was  noted  in  I  where  h(x)  = 
rectangular  wave  form. 
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We  observe  from  (36)  that  (30  does  depend  on  the  load  condition,  as  it  does  in  I.  Figure  3 
shows  (30  plotted  against  ?0  .  With  this  location  of  the  points  it  is  not  possible  to  determine  the 
precise  nature  of  the  functional  dependence.  However,  with  test  conditions  that  provide  well 
spaced  points  on  the  ?0  this  should  be  possible. 

It  follows  that  with  two  or  three  suitably  chosen  test  load  conditions  we  can  find  the 
average  value  of  c  and  obtain  a  graph  for  po  as  a  function  of  70  .  Then,  with  (9)  and  (20),  we 
can  predict  At;,  as  a  function  of  (Ayj.yj)  for  other  load  conditions  for  h(x)  =  sin  x,  e  =  1,  z  =  1/4, 
and  constant  or  variable  T,  and  also  find  tj  =  Atj  +  •  •  •  +  Atj_i  as  a  function  of  vj. 

We  now  take  up  a  number  of  interesting  points  not  previously  considered.  First,  using  the 
EHN  Data,  let  us  consider  the  effect  on  the  p,?0  estimates  when  e  takes  different  values  with 
the  same  5  =  .25  v.  We  only  consider  the  straight  portion  of  the  Pjsf0  .  Figure  4  shows  the 
results  for  e  =  ap/ae  -  1  equal  to  .25,  .50,  and  1.0.  The  straight  line  character  of  the  estimates 
still  obtains  in  each  case.  Reducing  e  causes  the  numerical  values  of  the  Pj?0  to  increase.  Our 
assumption  throughout  the  data  analysis  that  e  =  1  is  reasonable,  although  a  closer  connection 
with  how  the  form  of  s(t,  q)  evolves  with  time,  based  on  experimental  evidence,  is  certainly 
needed. 

Next,  let  us  hold  e  constant  and  change  6.  Figure  5  assumes  <=  =  1  and  8  is  changed.  As 
long  as  we  assume  5  =  zy,  we  see  that  the  straight  line  character  of  the  estimates  still  obtains. 
However,  if  5  =  constant,  the  straight  line  approximation  is  no  longer  satisfactory.  We  know 
that  the  length  of  the  high  stress  zone  in  the  crack  plane  increases  for  materials  that  work  harden. 
Thus,  the  choice  5  =  zy  is  reasonable  in  this  study.  Again,  experimental  evidence  on  how  s(t,  %) 
evolves  with  time  is  needed. 

Equations  (16),  (18),  and  (19)  contain  the  temperature  T  explicitly.  Let  us  consider  how  T 
influences  the  results  predicted  by  the  model.  We  employ  the  VIIG  Data  for  this  purpose. 


We  have  from  (3) 
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a  = 


(38) 


T0  =  307°K  for  the  data.  We  assume  s  =  1  and  z  =  1/4  and  employ  the  same  numerical  values 
for  the  other  quantities  as  previously  used  when  making  parameter  estimates  in  b). 

Let  Ti  =  312°K  and  T;  =  302°K.  We  find  using  (31)  and  (38)  that 

Tj  =3 12°K:  a  =  4.73250  x  KT25  ,  fys,,  '  =  25.5369  (1  -  .725106  yj) , 

T0  =  307°K:a=  1.15670  x  10"25  ,  (3;s0 '  =  25.9528  (1  -  .725106  yj) ,  (39) 

T2  =  302°K:  a  =  2.69906  x  10"26  ,  (3jT0  '  =  26.3826  (1  -  .725106  yj) . 

Notice  that  we  assume  the  coefficient  of  y  inside  parenthesis  and  also  X0  do  not  depend  on  T. 
We  use  (16)  and  Pj?0  from  (39)  to  predict  Atj  as  a  function  of  Avj  or  yj,  starting  at  crack  length 
16  mm  with  Aa;  =  1  mm  and  ending  at  crack  length  36  mm.  This  interval  16  mm.  -  36  mm.  is  in 
the  straight  portion  of  the  Pj?0  estimates.  We  could  have  started  at  9  mm.  if  we  were  willing  to 
employ  graphical  results  instead  of  (32). 

The  solid  curves  in  Figure  6  are  from  prediction  at  the  indicated  temperatures,  and  the 
solid  "dots"  are  from  the  data.  Two  points  are  worth  noting.  First,  as  we  expect,  the  model 
predictions  at  T  =  307°K  are  in  excellent  agreement  with  the  data.  Second,  the  predicted  Atj  at 
T=302°K  and  312°K  are  substantially  different  from  the  data  even  though  the  temperature 
difference  is  only  ±  5°K;  this  is  due  to  the  fact  that  while  the  temperature  difference  is  small  and 
the  change  in  the  p0?0  is  small,  small  differences  in  the  P0?0  produce  large  differences  in 
predicted  Atj  because  P0T0  appears  in  the  exponent.  We  do  not  have  data  with  which  to 
validate  either  prediction.  However,  the  experiments  of  Zhurkov  [3]  on  time  to  failure  under 
static  load  do  indicate  that  these  times  are  sensitive  to  small  temperature  changes.  The  results  in 
Figure  6  imply  that  unless  the  temperature  is  carefully  controlled  at  a  constant  value  or  is 
accurately  measured  throughout  the  experiment  the  data  will  be  of  little  value  for  model 
validation.  Further,  prediction  of  life  behavior  in  a  component  (regardless  of  which  model  is 
employed)  must  account  for  temperature  if  accurate  results  are  to  be  expected.  There  is  another 
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possible  aspect  of  temperature  behavior  that  we  must  consider  when  experiments  are  performed. 

We  speculated  in  Part  I  that  the  transient  behavior  observed  at  the  start  of  the  Pj?0 
estimates,  shown  for  example  in  Figure  2  might  be  due  solely  to  material  behavior.  However, 
examination  of  [5]  suggests  this  may  not  be  the  case. 

Considerable  care  was  taken  in  the  VHG  experiments  mean  of  by  a  suitable  load  shedding 
program  to  eliminate  the  effect  of  transient  material  behavior  due  to  the  rapid  starting  of  the 
crack.  Thus,  this  type  of  transient  material  behavior  is  not  likely  the  source  of  the  transient  Pj?0 
estimate  behavior  observed. 

We  note  that  in  the  VHG  experiments  desiccated  air  was  employed  and  the  room 
temperature  was  controlled  at  24°  C.  At  a  frequency  of  20  Hz,  we  know  the  specimen 
temperature  will  rise  above  24°  C  until  a  steady  value  is  reached,  and  we  assume  this  rise  is 
10°C.  giving  T  =  307°  K.  It  follows  that  in  the  VHG  Data  the  observed  transient  behavior  of  the 
P,J0  estimates  is  most  likely  due  to  the  transient  temperature  behavior,  this  also  is  most  likely 
the  case  for  the  IHN  Data  and  the  GD  Data.  However,  to  definitely  confirm  this  view  it  is 
necessary  to  measure  the  specimen  temperature  at  several  points,  particular  by  just  ahead  of  the 
crack  tip  in  future  experiments. 

Figure  7  is  concerned  with  the  effect  of  the  form  of  the  periodic  function  h(t)  on  the  Atj. 
To  study  this  effect,  we  again  employ  the  VHG  Data  with  T  =  307°K.  The  solid  curve  on  the 
left  is  for  h(x)  =  sin  r,  and  is  the  same  as  the  middle  curve  in  Figure  6.  We  have  no  experimental 
evidence  to  support  the  assumption  that  the  numerical  values  in  (32)  are  independent  of  the  wave 
form  in  h(t).  However,  let  us  assume  this  assumption  is  correct.  Let  us  assume  h(x)  is  given  by 
the  second  of  (13).  Then,  we  must  employ  (18),  and  predict  Atj  for  the  Ayj  from 

2Av-  ,  _ 

Atj  =  — ^-(ujAs  )e  Uj  ,  (40) 

CCZVj 

with  s  =  1,  z=  1/4,  and  the  PjS),  from  (32).  The  lower  solid  curve  in  Figure  7  is  predicted.  The 


difference  in  the  two  curves  is  substantial.  However,  since 


| sin  t |  >  I h(t) |  , 


where  h(r)  is  the  "sawtooth"  form,  the  difference  is  to  be  expected.  Again  we  have  no 
experimental  data  with  which  to  validate  the  prediction  shown  in  Figure  7. 

We  have  already  observed  that  small  changes  in  (30?0  due  to  small  temperature  changes 
produce  much  larger  changes  in  the  predicted  Atj  because  this  term  appears  in  the  exponent. 
Actually,  we  see  from  (40)  that  any  change  in  the  uj  given  by  the  first  of  (17),  even  by  small 
amounts,  produces  substantial  changes  in  the  predicted  Atj.  This  is  one  source  of  sensitivity  in 
the  model.  On  the  other  hand,  when  the  form  of  h(r)  is  changed  from  sin  x  to  "sawtooth,"  uj 
does  not  change;  rather  the  coefficient  of  e-u’  changes.  At  T  =  307°K,  we  find  for  the  VHG  data 
in  the  16  mm.  -  36  mm.  range  with  AUj  =  1  mm. 

(At.)Sawtooth  I - r 

— I - - - =  .797884 VujAs  ,  (41) 

(Atj)sutt  V  J 

and,  with  As  =  .401141,  and  Uj  -  52.5,  this  ratio  equals  3.66.  Thus,  the  predicted  Atj  also  are 
sensitive  to  the  reasonably  small  change  in  the  wave  form. 

Next,  let  us  observe  how  we  must  proceed  if  we  wish  to  predict  tj  with  different  ?Q  and  As. 
Consider  the  GD  Data.  It  is  clear  from  (37)  that  we  must  estimate  (30  and  c.  We  take  for  c  the 
average  value  .723957,  since  c  does  not  appear  to  depend  on  T0  and  As  >  0.  To  estimate  (30  for 
the  selected  Sg,  As,  and  ?0  =?0  +  As,  we  assume  that  it  suffices  to  approximate  the  (30  values  in 
(36)  by  the  curve 


6.714138 


.442001 

106 


.00401678 


For  the  selected  s0  ,  we  estimate  (30  from  this  equation  assuming  s0  is  either  in  or  not  to  far 
outside  the  interval  (6830,  10246)  given  in  (35).  We  then  employ  (16)  in  the  form 
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Av:  ,  ,, 

Atj  =  — j-  (2:uijAs  )h  e  1  , 
azvj 


(160 


and  Uj  from  (17)  and  (37)  with  e  =  1,  z  =  1/4  to  predict  the  tj  since  all  quantities  are  now  known. 

If  As  (or  As )  is  small,  the  approximation  (15)  will  not  be  satisfactory.  Then,  we  replace 
(16')  by 


Atj  = 


Avje  Uj 
azv; 


1-1 


pjraAs 


1  -  (1  ~  6  z)yj 


or,  more  generally, 


1  r-K  i-U-6  z)9i 

2tc  4  C 


(43) 


where  u-  =  Pj  ros0/{  1  -  €  z)y} }.  We  note  that  if  As  =  0,  i.e.  a  constant  stress,  (43)  because 

Av  e 

Atj*>  =  ^1-  ,  (44) 

azvj 

which  predicts  the  time  Atjs>  required  for  the  crack  to  increase  by  the  dimensionless  amount  Ayj 
under  a  static  far  field  tensile  stress  ?0.  Clearly,  the  Atjs)  in  (44)  will  be  much  larger  than  the  Atj 
in  (160  if  As  is  appreciably  greater  than  zero.  We  observe  that  the  second  factor  on  the  right  of 
(43)  represents  the  effect  of  As  in  reducing  Atjs)  to  Atj.  Since  we  know  (3j  depends  on  T0  and  As 
according  to  the  GD  Data,  and  since  we  have  no  data  for  the  static  load  case,  we  cannot  say  how 
the  Pj  in  (16  )  relate  to  the  Pj  in  (44). 

It  follows  from  the  material  presented  in  this  Section  that  the  data  analysis  and  study 
conducted  reveal  useful  information  about  the  model  and  its  predictive  capabilities,  and  the  FCG 
process.  It  also  is  clear  that  there  are  gaps  in  the  data  base  needed  to  validate  the  proposed 


deterministic  model. 
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i.  Discussion  of  Deterministic  Model 

We  consider  FCG  as  a  chemical  process  in  which  irreversible  atomic  bond  breakage  is 
responsible  for  the  growth  of  the  macro  observable  crack  length  and  for  which  rate  reaction 
theory  provides  the  basis  for  macro  modeling  starting  at  the  micro  level.  The  problem  we  face  is 
how  to  proceed  from  the  micro  level  to  the  macro  level  when  there  is  at  least  one  observable. 
Given  the  specimen  geometry,  a  periodic  tensile  load  F(t),  and  this  point  of  view,  let  us  consider 
the  bond  breaking  on  the  crack  plane  in  order  to  review  how  we  approach  this  problem. 

The  atomic  bond  breaking  process  on  the  crack  plane  is  complex.  An  elementary  picture 
of  this  process  is  as  follows:  In  the  immediate  region  to  the  right  of  the  crack  tip  at  a, 
irreversible  atomic  bond  breaking  due  to  the  external  load  is  the  dominant  process.  In  the  next 
region  on  the  right,  bonds  are  broken  by  the  external  load,  but  an  equivalent  number  of  new 
bonds  are  established  producing  plastic  deformation  as  the  dominate  process.  Only  the  thermal 
type  of  bond  breaking  and  healing  occurs  in  the  next  region  on  the  right;  this  is  the  region  of 
elastic  material  behavior.  We  model  this  elementary  picture  in  the  following  simple  manner. 

First,  we  assume  we  can  approximate  the  mean  bond  force  f  by  employing  the  average 
macro  concept  of  normal  stress  s(t,  g)  =  s(^)g(t),  where  g(t)  is  a  dimensionless  time  factor 
depending  on  F(t)  and  the  crack  plane  geometry,  and  s(^)  and  g(t)  are  bounded  functions. 
Second,  we  elect  to  divide  the  atomic  bond  breaking  process  into  a  region  just  to  the  right  of  the 
crack  tip  in  which  there  is  only  irreversible  bond  breakage,  and  a  region  to  the  right  of  the  first 
region  in  which  the  amount  of  'irreversible  bond  breakage  is  negligible,  disregarding  all  other 
types  of  bond  behavior.  Based  on  the  micro  result  (II)  from  reaction  rate  theory  for  the  mean 
rate 


AG*-Xf 

kT  kT  c  r  -m 
Kb  =  —  e  ,  f  =  s(t,  c)/v  . 

h 


(45) 


of  single  bond  breakage,  we  can  achieve  this  picture  of  the  process  by  taking  the  stress  factor 
s(g)  on  the  crack  plane  to  have  the  form  given  by  (5),  where  ap  >  ae  and  ap/Ge  =  ra  >  1.2.  The 
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expression  for  Kb  assumes  there  is  only  a  single  energy  barrier  that  must  be  crossed  for  atomic 
bond  breakage  and  there  is  no  bond  healing.  Clearly,  other  more  complex  Coemptions  can  be 
adopted  to  achieve  the  above  picture  of  the  bond  breaking  process;  further,  the  simplified  picture 
of  the  bond  breaking  process  can  be  enriched  in  various  ways. 

The  form  (4)  assumes  there  are  no  dynamical  effects  induced  by  the  testing  procedure,  and 
this  is  usually  the  case. 

The  macro  equations  (1)  and  (2)  are  based  only  on  our  picture  of  the  bond  breaking 
process  for  the  situauon  under  consideration.  On  combining  these  equations  with  (3),  (4)  and 
(5),  we  obtain  (6)  and  (7).  Equation  (7)  is  one  of  our  major  contributions,  and  is  a  different  form 
for  FCG  than  provided  by  other  approaches;  this  equation  is  interesting  for  several  reasons. 

First,  it  is  based  on  physical  principles  and  the  assumption  that  the  mean  bond  force  f  can 
be  approximated  by  employing  the  normal  stress  s(t,  q)  on  the  crack  plane. 

Second,  besides  the  dimensionless  crack  length  y,  the  quantities  5  and  rc  from  (5)  evolve 
with  time.  This  opens  up  the  possibility  that  y,  S,  and  ra  may  constitute  a  state  vector  if  we  can 
find  their  joint  evolutionary  equations.  This  possibility  has  not  previously  been  considered.  We 
do  not  have  such  joint  evolutionary  equations  at  this  time,  however.  Thus  this  interesting 
possibility  cannot  be  considered  further. 

Third,  (7)  contains  explicitly  the  temperature  T,  and  the  load  F(t)  through  the  far  field 
tensile  stress  s0(t).  Thus,  (7)  contains  quantities  of  direct  engineering  interest  that  help  define 
the  environment. 

Fourth,  (7)  contains  explicitly  the  material  parameters  a  (  or  AG+)  and  [3  as  defined  by  (3). 
AG~  is  a  free  energy;  it  is  a  thermodynamic  quantity  whose  value  either  can  be  estimated  once 
for  all  for  the  material  from  tabulated  data  on  the  heat  of  sublimation  (see  Zhurkov  [3]  and 
Krausz-Evring  [8])  or  it  can  be  estimated  from  simple  static  tests  (see  Zhurkov  [3]).  On  the 
other  hand,  the  growth  control  material  parameter  p,  which  we  find  depends  on  y,  defines  the 
strength  of  the  material;  it  depends  on  how  the  material  is  formed  into  components  and  must  be 
estimated  by  means  of  appropriate  test  data.  Thus,  (7)  contains  material  parameters  of  direct 
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engineering  interest. 

Fifth,  the  alternative  form  of  s(c)  given  in  (17bl)  admits  the  possibility  that  the  highest 
stress  is  at  the  crack  tip.  This  means  that  even  if  Act  is  small  the  maximum  rate  of  bond 
breakage  will  occur  at  the  crack  tip,  which  is  physically  appealing.  This  point  is  discussed  in 
Part  I  (see  eq.  (19bl)),  where  we  note  that  if  Ac/ce  is  small,  the  form  i  f  (7)  remains  unchanged. 
Thus,  (7)  still  applies  in  this  more  general  and  physically  appealing  situation. 

Sixth,  inspection  of  the  integral  in  (12)  suggests  that  there  might  be  different  h(x)  that  will 
give  the  same  value  to  the  integral;  i.e.  different  wave  forms  will  produce  the  same  damage 
accumulation  process  if 

_Lr-*e*h(T)dx 

2a:  ■'o 


is  the  same  for  different  h(x)  and  all  x  >  0.  Given  the  conditions  assumed  for  h(x)  plus  the 
continuity  conditions  imposed  by  the  fact  that  it  is  a  physical  load,  it  can  be  shown  this  is  not 

possible,  if,  in  addition,  h(t)  is  monotone  increasing  (or  decreasing)  on  the  interval  y)  as 

is  the  case  when  h(x)  is  given  by  (13).  If  this  monoticity  condition  is  not  satisfied,  however, 
different  h(x)  will  produce  the  same  value  for  the  integral;  while  this  situation  has  interesting 
implications,  we  will  not  discuss  it  here. 

Finally,  if  there  is  a  chemical  reaction  involved,  such  as  corrosion,  it  will  enter  (7)  through 
the  growth  control  parameter  {3.  Now  let  us  turn  to  how  (7)  is  employed  in  conjunction  with 
FCG  data. 

The  integration  of  (7)  requires  that  we  assume  a  specific  form  for  s0(t);  we  take  sQ(t)  to  be 
tensile  and  periodic  in  t  as  specified  in  (8).  Further,  since  we  do  not  have  joint  evolutionary 
equations  for  y,  5,  and  r0  at  this  time,  we  assume  the  constraint  equations  (11)  for  5  and  rc.  The 
combination  of  (7),  (8),  and  (11)  produces 


28 


dv 

— -  ctzve 
dt 


{3yrg(s0  Chicot)) 
1 -(1-<S  z)v 


(46) 


since  if  ra  >  1.2,  the  second  term  on  the  right  of  (7)  can  be  ignored.  Even  with  these  choices,  the 
exact  integration  of  (46)  is  difficult.  However,  the  use  of  different  forms  for  the  periodic 
element  h(t)  in  (8),  for  example  (13),  makes  possible  the  approximate  integration  of  (46)  in  the 
difference  forms  (16),  (18),  and  (19)  for  the  deterministic  development  of  FCG  in  discretized  t 
and  v;  this  step  is  reasonable  since  y  only  changes  by  small  amounts  for  all  Atj  involved. 

The  selection  of  the  constraint  equations  is  based  on  what  we  know  about  the  general 
behavior  of  FCG  in  ductile  polycrystaline  metals,  as  already  noted,  and  simplicity. 
Experimental  evidence  must  ultimately  be  employed  to  determine  if  the  choice  is  consistent  with 
the  facts.  Further,  if  (11)  is  a  reasonable  choice,  it  is  then  necessary  to  determine  on  the  same 
basis  if  our  choice  £  =  1  and  z  =  1/4  is  reasonable.  Such  experiments  would  be  different  from 
those  usually  conducted  in  FCG.  In  addition,  these  experiments  would  have  to  measure  the 
temperature  T  as  we  have  found  it  to  be  important  parameter. 

We  now  wish  to  comment  on  two  other  models  of  FCG  that  originate  in  reaction  rate 
theory.  Yokobori  and  co-workers  [10]  have  since  the  mid  1950’s  applied  reaction  rate  theory  to 
fracture,  fatigue,  and  FCG,  producing  many  interesting  results;  we  note  in  particular  [11]. 
Another  effort  in  this  direction  for  FCG  is  given  in  [12].  Each  introduces  crack  length  and  load 
in  a  different  manner  than  we  do.  In  [12],  for  example,  fracture  mechanics  based  quantities  are 
introduced  in  place  of  the  Xf-term  in  the  exponent  of  (II)  or  (45)  which  is  physical  questionable. 
Each  contains  one  material  parameter  equivalent  to  our  a.  Each  contains  a  second  material 
parameter  similar  to  our  (3;  however,  neither  group  investigated  to  determine  if  this  parameter 
depends  on  crack  length.  We  have  found  in  our  formulation  that  [3  depends  on  y  and  this 
dependence  is  essential  for  predictive  accuracy.  It  therefore  appears  that  our  method  for 
introducing  crack  length  and  load  (stress)  has  advantageous  not  possessed  by  these  two  models. 

We  assume  the  periodic  form  (8)  for  s0(t)  because  this  is  the  case  for  the  available  data. 
However,  it  is  clear  that  in  (7)  there  is  no  mathematical  restriction  on  sQ(t)  beyond  being  a 
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positive  bounded  cornin'  jus  function,  oscillatory  or  not.  For  example,  s0(t)  can  be  an 
oscillatory  non-periodic  function  such  as  provided  by  a  sample  function  from  spectrum  loading. 
But,  we  must  then  confront  the  possibility  that  hi-lo  vs.  lo-hi  loads  may  produce  different 
amounts  of  physical  damage  accumulation  than  predicted  mathematically  by  (7)  due  to  the  load 
history  dependence  in  the  material  [9].  Put  another  way,  the  growth  control  parameter  p  may 
depend  on  load  order.  Without  experimental  evidence  to  lend  support,  predictions  of  FCG  based 
on  (7)  for  any  such  s0(t)  must  be  viewed  with  caution. 

The  model  consider  so  far  is  deterministic.  The  data  analysis  employs  the  pair  of 
quantities  Avj  and  Atj.  Ayj  is  specified  by  the  choice  of  the  division  aj  of  the  crack  length 
employed.  The  time  tj  to  reach  aj  depends  on  the  sample  selected  in  a  test,  and  these  times  vary 
in  a  random  manner.  Let  the  random  variable  Tj  denote  this  random  time.  We  assume  for  a 
deterministic  model  that  Atj  =E(ATj};  i.e.  the  deterministic  model  only  applies  to  mean  time 
behavior.  Even  with  the  sample  sizes  employed,  the  E{ATj}  do  not  decrease  in  a  smooth 
manner.  Thus,  there  is  a  fog  imposed  by  the  finite  sample  size  and  the  experimental  technique 
through  which  we  must  examine  results  obtained  in  the  data  analysis,  the  smaller  the  sample  size 
and  the  poorer  the  experimental  technique  the  denser  the  fog. 

Equation  (16)  is  used  to  estimate  the  Pj70.  In  data  sets  b),  c),  and  d),  the  estimates 
approximate  a  straight  line  except  for  the  first  few  initial  points;  thus  transient  behavior  is 
initially  present  before  regular  behavior  is  observed.  We  have  already  indicated  that  this 
transient  behavior  is  most  likely  due  to  the  small  gradual  rise  in  specimen  temperature  resulting 
from  internal  material  hysteresis  before  a  constant  temperature  slightly  above  the  ambient  air 
temperature  is  reached. 

We  also  observe  that  the  straight  line  character  of  the  pjT0  estimates  persists  to  the  end  of 
the  crack  length  generated  even  though  in  some  cases  the  final  crack  length  is  close  to  the 
specimen  boundary.  Thus,  the  effect  of  the  boundary  does  not  seem  to  be  significant. 

In  all  events,  we  find  for  each  of  the  six  data  sets  considered  the  straight  line  form 
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1*0  Po*0  (1  V.Vj) 

(47) 

v*o  Po*o  (1  cy) 

(48) 

where  (30  =  A0/ukT.  With  crack  length  as  the  only  macro  observable,  we  find  there  are  three 
material  parameters  present,  namely  a  or  AG+,  j30  or  X^/x),  and  c.  Zhurkov  and  Krausz-Eyring, 
employ  only  two  constant  material  parameters  in  their  model  of  time  to  failure  under  constant 
load,  and  static  load  data  support  their  parameter  assumptions.  In  some  fatigue  life  models  (no 
observable)  [14,15],  two  constant  material  parameters  also  are  employed.  Further,  the  models  of 
FCG  presented  in  [11,12],  where  there  is  an  observable,  also  assume  two  constant  material 
parameters;  however,  these  investigators  did  not  perform  the  same  type  of  evolutionary  (yj,  tp 
data  analysis  as  we  employ,  and  thus  it  is  not  possible  to  determine  if  their  parameters  depend  on 
v.  In  the  model  of  FCG  presented  here,  it  is  essential  for  evolutionary  predictive  accuracy  that 
there  be  three  material  parameters.  To  see  this,  let  us  write  (e  =  1,  z  =  1/4,  h(x)  =  sin  x) 


-PjSp '  _  2(30s0,(l  -cyj) 

1  -  .75  vj  1  -  .75  yj 


(49) 


which  is  the  only  quantity  involving  the  two  material  parameters  p0  and  c.  If  does  not  depend 
of  y,  then 


2PoS0' 
l  “  -75  vj  ' 


(50) 


Let  us  calculate  Atj  from  (16)  using  these  equations  and  the  VHG  Data.  At  yj  =  .216535 
(aj  =  16.5  mm.),  let  us  fix  the  pQ?0  in  (50)  so  that  the  (Atp2  obtained  using  (16)  is  the  same  as 
the  (Atj)L  obtained  using  (49)  in  (16).  Then,  at  yj  =  .464878  (aj  =  35.5  mm.),  we  find  that  the 


ratio 
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(At;)!  , 

— =  6.135490  x  1  O'6  .  (51) 

( Atj  )2 

It  follows  that  if  we  neglect  the  dependence  of  pjS0  on  vj,  as  is  done  in  (50),  we  will  predict  tj 
that  are  much  smaller  than  observed.  Thus,  we  see  that  if  only  time  to  failure  is  of  interest,  then 
two  constant  material  parameters  appear  to  suffice,  but  if  evolution  of  crack  growth  is  of  interest 
three  constant  material  parameter  must  be  used. 

We  know  from  Part  I  that  a  decrease  in  X  means  an  increase  in  material  strength.  Since 
from  (37)  Xy  =  XQ(l  -  cy),  it  is  clear  that  the  material  is  increasing  in  strength  at  a  constant  rate 
with  respect  to  y,  which  implies  the  materials  employed  show  work  hardening.  This  implication 
is  consistent  with  what  is  known  about  the  behavior  of  such  materials.  What  do  we  currently 
know  about  what  P0  or  X0  and  c  are  dependent  upon? 

First,  let  us  note  some  information  provided  by  other  investigators.  The  data  in  [3]  suggest 
that  a0  and  c  do  not  depend  on  T.  Further,  [16]  suggests  that  X0  and  c  depend  on  how  the 
material  is  formed  into  components.  We  shall  assume  both  hold.  However,  direct  experimental 
evidence  is  needed  to  verify  these  assumptions. 

P0  is  related  to  a.0  through  the  second  of  (38);  i.e.  p0  depends  on  T  although  we  are 
assuming  X0  does  not.  We  shall  confine  attention  to  PQ  in  what  follows. 

The  analysis  of  the  GD  Data  shows,  by  means  of  (36)  and  Figure  4,  that  PQ  depend  on  ?0 
or  s0  and  As,  assuming  s  =  1,  z  =  1/4,  1  hit)  =  sin  x.  Figures  5  and  6  and  their  ms. 

approximate  straight  lines  reveal  that  P0  also  depends  on  e  and  z  when  h(t)  =  sin  i.  Thus,  given 
the  data  analyzed  and  the  material  employed,  we  know  that  P0  depends  on  "s0,  As,  T,  e,  and  z 
when  h(x)  =  sin  t.  However,  we  do  not  have  data  with  which  to  determine  if  PQ  depends  on  the 
periodic  wave  form  h(t)  for  fixed T0,  As,  e ,  z,  and  T. 

t 

Actually,  in  (42)  we  assume  P0  depends  on  just  s0  .  This  is  a  reasonable  assumption  since 

the  dependence  of  P0  on  As',  is  relatively  weak  in  comparison  to  its  dependence  on  T0  for  As 

/ 

not  small.  Even  so,  the  values  of  s0  employed  in  (42),  and  by  implication  the  values  of  s0  and 
As,  should  not  be  too  far  outside  the  range  of  values  employed  in  the  GD  Data.  Next,  consider 
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what  we  have  learned  about  c. 

The  adjacent  table  has  been  complied  from  results  given  in  the  previous  section. 


Table  1 


Material 

s0 

Po 

c 

Al.  alloy  2024-T3 

8767 

.00296039 

.725106 

Al.  alloy  2024-T3 

10820 

.00246603 

.723869 

Al.  alloy  7075-T6 

10246 

.00260708 

.722287 

Al.  alloy  7075-T6 

10004 

.00269365 

.727481 

Al.  alloy  7075-T6 

6830 

.00388267 

.722104 

The  remarkable  point  suggested  by  this  table  is  that  irrespective  of  the  aluminum  alloy, 
specimen  thickness,  and  the  ?0  considered,  c  remains  essentially  the  same  with  an  average  value 
.724051  and  a  maximum  deviation  of  .47%  for  €  =  1,  Z=  1/4,  different  temperatures,  and 
h(t)  =  sin  t.  More  data  are  needed  to  investigate  this  point. 

We  know  from  Figures  5  and  6  and  the  associated  m.s.  approximate  straight  lines  that  c 
depends  on  s  and  z.  Other  conditions  being  the  same,  we  have  no  data  on  how  c  might  vary 
with  h(x).  Thus,  based  on  the  specimen  geometry  and  material  considered,  c  depends  on  e  and  z 
and  possibly  h(x),  but  does  not  depend  on  I0  and  T.  These  observations  have  implication 
concerning  the  amount  of  testing  required  for  accurate  prediction. 

Accurate  evolutionary  prediction  under  different  load  conditions  and  temperature  requires 
that  we  have  accurate  estimates  of  P0  and  c  as  a  function  of  J0  for  a  given  component  geometry, 
T,  € ,  z,  and  wave  form  h(t),  and  assuming  no  transient  material  behavior.  If  we  have  three 
reasonably  high  replication  sets  of  FCG  tests  with  the  sQ  substantially  different  but  with  the 
other  quantities  fixed,  we  can  estimate  P0  or  X^  and  c  for  each  test  s0.  We  plot  p0  or  X^  and  c 
against  ?0  and  put  a  smooth  curve  through  these  points.  At  any  other  1Q  (not  too  far  out  of  the 
range  of  the  test  values),  we  can  obtain  appropriate  values  for  P0  or  /i0  and  c.  Then,  employing 
(16),  we  can  predict  mean  t,  for  this  70  and  T  if  different  from  the  test  values.  While  these  tests 
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are  not  likely  required  for  obtaining  c,  we  know  they  are  for  estimating  p0  or  X0,  assuming  Xq 
and  c  do  not  depend  on  T. 

Equation  (44)  merits  a  few  comments.  First,  as  must  be  the  case,  we  obtain  this  equation 
when  using  the  other  two  wave  forms  provided  we  correct  the  right  based  side  of  (18)  and  of 
(19)  so  that  As  >  0  can  be  small.  Second,  let  us  rewrite  (44)  as 


Av 


At, 


ra3J.'(l-cyj) 
j  —  1— (1 — e  z)y, 

-  =  cczvie  7) 


(440 


We  can  now  compare  (440  with  (3)  and  (10)  in  [13],  on  taking  into  account  differences  in 
notation.  We  observe  in  this  static  load  case  the  results  differ  because  we  are  following  the 
evolution  of  an  observable  with  a  nonuniform  stress  distribution  and  they  do  not  do  this. 

Let  us  rewrite  (19)  as 

=  H>-“*  ,  (191 

-  azv: 

-  j 


which  is  the  same  form  as  (44),  except  for  a  factor  two,  even  though  in  (19)  we  have  a  periodic 

_ _  * 

load  and  As  is  not  small.  It  is  clear  that  when  h(x)  has  a  rectangular  wave  form  and  sQ  =  s0  ,  the 
static  load  produces  the  same  evolutionary  result  in  half  the  time  provided  the  pj  remain  the 
same.  This  possibility  was  noted  by  Regel  and  Leksovsky  [15]  although  their  approach  to  the 
fatigue  problem  is  very  different  from  ours  but  similar  to  that  employed  in  [14], 

Finally,  we  note  that  it  is  possible  to  use  the  deterministic  model  for  accelerated  testing  in 
useful  ways.  One  way  to  accelerate  a  test  is  to  either  increase  ?0,  or  As,  or  both  by  constant 
amounts  over  the  condition  at  which  a  prediction  is  required.  We  next  estimate  P0  and  c  at  the 
accelerated  test  condition.  If  we  assume  these  estimates  will  apply  at  the  10  and  As  for  which 
prediction  is  required,  then  we  can  employ  them  to  produce  the  predicted  behavior.  However, 
we  know  that  P0  depends  on  s0  and  As;  i.e.  the  load  interaction  will  be  different  at  the  two  load 
conditions.  The  same  problem  was  noted  in  [20].  Hence,  the  accuracy  of  the  predicted  behavior 
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is  open  to  question  when  T0  and  As  are  changed  to  accelerate  the  test.  However,  since  a  and  [3 
depend  on  T  in  a  known  manner,  another  way  to  accelerate  a  test  is  to  increase  T.  The  second 
way  has  some  advantages. 

We  already  know  from  Figure  7  and  from  [3,  16]  that  a  small  change  in  T  produces  a 
significant  change  in  mean  life  behavior.  Suppose  we  want  to  predict  mean  life  behavior  at 

o  _  o 

T0  =  302  K  with  specified  So,  As,  <=  ,  z,  and  h(t)  =  sin  t.  Let  us  test  at  T2  =  312  K  with  the  same 
specified  quantities.  Figure  7  indicates  that  the  time  of  specimen  testing  is  reduced  by  a  factor  of 

o 

about  3.  If  we  test  instead  at  T3  =  322  K,  the  time  the  specimens  are  under  load  is  reduced  by  a 
factor  of  about  10.  This  method  of  accelerating  testing  has  the  advantage  of  eliminating 
questions  that  may  arise  concerning  the  effect  on  mean  life  behavior  due  to  load  interaction 
when  changing  and  As  with  T  held  constant.  The  second  wav  of  accelerating  a  test  thus  has 
advantages  over  the  first  way. 


5.  Probabilistic  Model 


The  model  discussed  so  far  is  deterministic.  However,  the  FCG  phenomena  is  statistical. 
Even  with  controlled  temperature  T,  load  F(t),  and  humidity,  as  in  the  VHG-experiments, 
variability  among  the  results  for  different  test  specimens  is  significant;  the  variability  in  this  case 
is  primarily  due  to  variability  in  material  behavior.  A  second  and  different  source  of  variability 
in  FCG  is  the  variability  always  present  in  the  load,  temperature,  and  environment;  i.e.  in 
experimental  technique.  These  are  the  sources  of  the  statistical  nature  of  the  FCG  phenomenon. 
Inspection  of  the  sample  function  graphs  of  the  data  sets  of  crack  growth  in  [17]  indicates  that 
the  VHG  Data  has  the  least  scatter.  This  implies  that  the  noise  due  to  experimental  technique  is 
least  for  these  data.  Thus,  we  expect  the  least  scatter  in  the  moment  estimates  and  edFs 
(empirical  distribution  functions)  of  the  ATj  and  Qj  (see  below)  in  the  VHG  Data,  and  this  is  the 
case.  In  this  paper,  we  employ  an  elementary  approach  in  our  first  attempt  to  model  material 
variability. 

Consider  (16): 


Av-  V/^Tt  / 

At.  =  E  { AT . )  =  ~  AiUj  As'  e“Uj  ,  (52) 

J  J  CCZV;  J 

'  J 


where  we  specifically  note  h(t)  =  sin  t.  As  is  not  small,  and 


WT 


Pj(l  +  6)Sp 

1  -  (1  -  6  z)Vj  ’ 


We  assume  T0,  As,  T,  and  the  environment  are  controlled  to  be  constant.  We  have  already 
assumed  in  these  equations  that  the  mean  bond  force  f  is  specified  by  the  deterministic 
expression  f  =  s(t,  q)/u.  Thus,  the  only  sources  of  material  variability  lie  in  AG+  and  PjS),  . 
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AG~  is  a  free  energy,  for  such  a  thermodynamic  quantity,  the  variability  is  small  so  that  in 
this  initial  stage  of  modeling  variability  we  will  neglect  it  even  though  a  is  sensitive  to  changes 
in  AG+.  Thus,  (3,?0  is  the  only  source  of  variability  that  we  will  consider  at  this  time. 

Let 


zj  =  a/ujAs  e  ; 


it  follows  from  (52)  that 


Atj=E{ATJ}  = 


AyjVijr 

ocZyj 


Next,  let 


Ayi 

ATj  =  ■  —  Zj  , 


OCZV; 


where  =  (1  +  Pj)zj  and  E  { Pj }  =  0.  Then,  from  (54), 

Av;  ViriT 

E(  ATj }  =  -  -J-  —  Zj, 

OZV; 

*  J 

which  is  (52). 

Now  let  Uj  and  Qj  be  two  additional  rvs  with 

Uj  =  Uj (1  +  Qj) , 


and  let 


(53) 


(54) 


Zj  =  a/UjAs  e  Uj  =  a/uj(1  +  Qj)As  e  U)(1  + 


(55) 


We  now  have 
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E 


/U;As  e  Uj 


uj(l+Qj)As  e 


-Uid  ■ 


and  hence  it  follows  that 


E  (1+Qj)1/2  e 


-uiQi  U 


=  1 . 


(56) 


The  quantity  uj  is  of  the  order  of  50  for  Al..  This  means  that  the  Qj  are  small.  Then  with 

log  Zj/Zj  =  i  log  (1  +  Qj)  -  UjQj 


it  follows  that 


Zj/zj 


-(Uj-yKJj 

e 


log  Zj/zj 


(57) 


At  the  sample  function  level 


qj.n- 


lOg  Zj,n/Zj 


(58) 


where 


zj.n/zj  - 


^^j.n 

eTat -} 


(59) 
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Further, 


p  Qj^qj  =p 


log  Z;/z; 

- i — i.  <  n- 


sqj 


Uj - 


=  p 


Zj/zj  >  e 


.  1  > 


K  or 


FQj(qJ)=l-FZj/2i 


-frrytoi 


(60) 


This  method  of  introducing  the  Qj  avoids  addressing  the  joint  problems  of  how  to  assign 
variability  at  the  micro  material  level,  and  how  to  convert  such  a  variability  into  macro 
stochastic  FCG  equations  in  a  physically  satisfactory  manner.  We  defer  addressing  these  two 
problems  at  present,  and  go  on  to  determine  what  a  data  analysis  says  about  the  statistical 
properties  of  the  Qj.  We  first  consider  the  VHG  Data  which  has  the  highest  replication  number 
(sixtv-eight)  and  the  variability  at  the  sample  level  is  small.  We  assume  e  =  1,  z=  1/4,  and 
h(x)  =  sin  t. 

First,  consider  the  four  crack  length  intervals  of  different  lengths  centered  as  closely  as  the 
data  permits  on  the  crack  length  26  mm: 

(25.8  mm.,  26.2  mm.),  (25.6  mm.,  26.6  mm.),  (25  mm.,  27  mm.),  and  (24.6  mm.,  27.6  mm.) ; 

(61) 

The  crack  lengths  of  these  intervals  are  respectively.  .4  mm.,  1  mm.,  2  mm.,  and  3  mm..  We 
find  from  the  data  in  each  case  the  sample  times  Atjn  required  to  cross  these  intervals,  where 
n  =  1,  •  •  ■  ,  68  and  j  =  1,  •  ■  • ,  4  refers  to  the  above  intervals.  We  next  estimate  the  E{ATj }  and 
compute  the  sample  qj  n  from  (58).  The  estimated  moments  of  the  Qj  are  given  in 
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Table  2 


j 

rriQ 

°Q 

mQ/aQ 

Vi.j 

1 

256.265  x  10~6 

9.324534  x  10"6 

.083922 

-.45514 

2 

168.294  x  10"6 

6.099563  x  10“6 

.068143 

-.355644 

j 

123.985  x  10~6 

4.521016  x  10"6 

.058313 

-.631304 

4 

98.6324  x  10"6 

3.535632  x  10"6 

.052455 

-.909978 

The  smoothed  edf’s  of  the  standardized  variable  qj,n  =  (qj,n  ~  mj)/Cj  are  plotted  in  Figure  8. 
These  smoothed  edf’s  am  obtained  from  the  smoothed  edf’s  of  the  ATj  employing  (60).  The  qj>n 

-j  -j 

range  from  approximately  -8  x  10  to  7  x  10  . 

The  means  are  small  and  decrease  as  the  crack  length  interval  increases.  The  variances 
decrease  as  the  crack  length  interval  increases.  These  results,  particularly  the  last,  suggest  that 
the  Qj  are  the  result  of  some  type  of  smoothing  operation  over  increasing  crack  lengths  at  the 
micro  level. 

The  skewness  coefficients  Yi.j  are  all  negative  with  average;  -.588110  and  not  close  to 
zero.  This  results  indicates  that  the  Qj  are  not  normal  random  variables. 

That  the  smoothed  edf’s  are  close  is  not  surprising  in  view  of  the  fact  that  this  is  the  case 
for  the  smoothed  edf’s  of  the  Atj  (see  [16]). 

Next  consider  the  four  crack  length  intervals 

(21  mm.,  22  mm.),  (26  ntm.,  27  mm.),  (31  mm.,  32  mm.),  and  (36  mm.,  37  mm.) , 


each  having  Aaj  =  1  mm..  The  estimated  moments  of  these  Qj  are  given  in 
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Table  3 


j 

mQ 

Sq 

m  q/Cq 

Yi.j 

1 

189.706  x 

10-6 

6.200268  x  10"6 

.076186 

-.07860 

2 

89.9706  x 

10~6 

5.088137  x  10"6 

.038556 

-.372504 

'i 

j 

190.265  x 

10"6 

7.719240  x  10"6 

.068481 

-.547555 

4 

181.074  x 

10"6 

9.981985  x  10~6 

.057312 

-.284855 

The  correlation  coefficients  pjk  ire  estimated  and  given  in 

Table  4 


j.  j  +  1 

j.j  +  2 

j.  j  +  3 

.154233 

.132162 

.370998 

.173791 

.240113 

.470985 

The  smoothed  edf’s  of  the  standardized  qj  n  are  shown  in  Figure  9,  with  the  same  approximate 
range  for  the  qj  n  as  in  the  previous  case. 

The  estimated  means,  variances,  and  skewness  coefficients  fluctuate  about  a  constant 
values;  the  averages  of  these  estimates  are  162.004  x  10-6,  7.247  408  x  10-6,  and  -.320879, 
respectively.  Again  the  nonzero  skewness  coefficients,  which  are  negative,  indicates  that  these 
Qj  are  also  not  normal  random  variables.  The  correlation  coefficients  are  all  positive;  thus,  the 
Qj  are  statistically  dependent.  The  smoothed  edfs  of  the  standardized  random  variables  are  also 
close  to  one  another  and  the  set  of  edfs  in  Figure  8  is  close  to  the  set  in  Figure  9.  These  results 
are  for  one  ?0  and  As  with  different  crack  length  intervals. 

Next,  let  us  consider  the  GD  Data  where  there  are  three  load  conditions  (s0.  As),  and  there 
are  sixty  replications  for  each  condition.  We  consider  three  crack  length  intervals: 
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(16  mm.  -  17  mm.),  (17  mm.  -  18  mm.),  (16  mm  -  18  mm.) 

The  results  presented  here  are  for  different  loads  but  for  each  of  the  indicated  crack  length 
intervals. 

Figures  10,  11,  12  display  the  edf  s  of  the  standardized  Q’s  for  the  three  load  conditions  at 
the  indicated  crack  length  intervals.  These  edf’s  are  obtained  from  the  edf’s  of  the  ATj  by 
employing  (55)  and  (60).  We  see  from  the  figures  that  the  edf’s  of  the  standardized  Qj  are  close 
to  one  another  for  different  s"0  and  As.  However,  we  see  from  Tables  5,  6,  and  7  that  some  of  the 
y’s  are  positive;  this  is  likely  due  either  to  sample  size  variability  or  to  variability  in  the 
experimental  program  or  to  both,  and  occurs  when  the  positive  ’s  of  the  ATj  are  much  less 
than  one. 

Tables  5,  6,  and  7  also  give  estimated  moment  results  for  Q’s  at  the  same  crack  length 
intervals.  We  observe  that  the  moments  of  the  Q’s  not  only  depend  on  the  crack  length  interval 
but  also  on  the  load.  The  Q’s  for  the  intervals  (16  mm.,  17  mm.)  and  (17  mm.,  18  mm.)  are 
correlated  as  is  the  case  for  all  non-overlapping  intervals,  which  is  not  surprising  since  we 
already  know  the  ATj  on  non-overlapping  intervals  are  positively  correlated.  We  observe  from 
these  tables  that  the  mean  mq,  the  variance  a§,  and  the  signal  to  noise  ratio  mq/oq  increase  in 
all  cases  as  ?0  decreases.  However,  since  the  T0  are  not  well  spaced,  we  cannot  form  a  firm 
idea  of  how  these  quantities  depend  on  70  ;  this  point  must  be  taken  up  at  a  later  time  when  data 
are  available  at  more  than  three  points  or  when  there  is  better  spacing  in  the  J0  values. 

Table  5 


Q,  (16  mm.,  17  mm.) 


mQ 

aQ 

mq/CQ 

Yi 

$0 

SI 

.292595  x  1(T3 

.108816  x  10-4 

.0886993 

-.294657 

10246 

S2 

.398630  x  10"3 

.146009  x  10-4 

.104323 

-.430016 

10004 

S3 

.984548  x  10“3 

.374576  x  10~* 

.160867 

-.0755493 

6830 
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Table  6 


Q,  (17  mm.. 

18  mm.) 

mQ 

ctQ 

mq/CQ 

Yi 

s</ 

SI 

.222031  x  10“3 

.851427  x  10"5 

.0760923 

.210108 

10246 

S2 

.363674  x  10"3 

.135120  x  10-4 

.0989356 

-.249759 

10004 

S3 

.942224  x  10"3 

.383306  x  10-4 

.156321 

.00992309 

6830 

Table  7 

Q,  (16  mm.  - 

18  mm.) 

ittq 

aq 

TDq/Oq 

Yi 

s</ 

SI 

.229298  x  10~3 

.876723  x  10"5 

.0776100 

.111876 

10246 

S2 

.356154  x  10*3 

.130594  x  10~* 

.098554 

-.426007 

10004 

S3 

.880428  x  10~3 

.340063  x  lCT4 

.150975 

.0728710 

6830 
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6.  Discussion  of  Probabilistic  Model 

We  have  introduced  randomness  into  the  macro  deterministic  difference  form  (16)  of  the 
FCG  process  by  the  expedient  of  assuming  the  growth  control  parameter  is  muluplied  by  a 
random  factor.  We  have  already  commented  on  the  fact  that  this  method  of  introducing  the  ever 
present  variability  may  have  short-comings  when  the  physics  of  the  random  micro  material 
behavior  is  the  starting  point  for  introducing  randomness.  The  behavior  of  mean  and  variance  of 
the  Qj,  when  the  crack  length  intervals  (61)  are  employed  in  the  VHG  Data,  suggests  that  these 
Qj  are  the  result  of  the  some  smoothing  operation  at  the  micro  level  over  the  intervals,  which 
indicates  the  micro  starting  point  merits  attention.  Further  investigation  on  this  approach  is 
required. 

The  closeness  of  the  smoothed  edf  s  of  standardized  Q’s  in  Figures  8  and  9  for  different 
crack  length  intervals  but  the  same  ?0  and  As  is  similar  to  what  was  reported  earlier  [17,  18]  for 
the  ANj  (ATj  in  present  paper).  What  is  of  interest  is  that  the  edf  s  of  the  standardized  Qj  are 
close  to  one  another  for  the  same  crack  length  interval  but  with  different  ?0  and  As,  as  is  clear 
from  Figures  10,  1 1  and  12  obtained  from  the  GD  Data;  further,  if  one  compares  the  results  in 
this  set  of  figures,  we  observe  that  these  edf  s  are  close  to  one  another  even  though  the  crack 
length  intervals  are  different.  All  of  this  suggests  that  the  standardized  Q;  appear  to  be  the  same, 
having  a  negative  skewness,  in  spite  of  the  fact  we  are  viewing  them  through  the  fog  imposed  by 
the  finite  sample  size  and  experimental  variability. 

The  correlation  coefficients  in  Table  4  are  not  large  but  are  positive.  Since  the  intervals  in 
(62)  are  well  separated,  this  results  is  not  surprising,  especially  in  view  of  the  results  on  the 
correlation  coefficients  of  the  ANj  given  in  Table  3  of  [17]  and  first  reported  and  modeled  by  us 
in  [18],  Thus,  the  statistical  dependence  of  the  Qj  is  to  be  expected.  This  means,  of  course,  that 
the  Qj  form  a  family  of  dependent  random  variables  indexed  on  the  interval  Ayj;  i.e.  they 
constitute  a  random  function. 

The  non  normal  character  of  the  Qj  can  be  explained  as  follows:  in  virture  of  the  fact  that 
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(1  +  €  )(1  +  Qj)PjS0 

U  +Qj)Uj  = - — — , 

l-(l-e  z)vj 

we  observe  that  the  (1  +  Qj)(3jJ0  represent  the  random  nature  of  the  dimensionless  work  done  in 
breaking  a  bond.  Since  this  work  depends  on  the  variability  inherent  in  the  bond  force  f  and  in 
the  distance  A.  a  bond  is  extended  in  breaking,  the  variability  in  the  work  done  in  breaking  a  bond 
is  dependent  on  the  joint  distribution  of  these  two  random  elements.  Thus,  the  non  normal 
character  of  the  Qj  is  not  unexpected.  We  also  note  that  the  negative  skewness  of  the  Qj  or 
1  +  Qj  suggests  the  possibility  that  the  origin  of  these  rv’s  might  be  associated  with  a  "weakest 
link"  concept.  In  the  present  context,  this  suggest  that  the  underlying  bond  breaking  process 
might  be  associated  with  (governed  by)  a  mechanism  based  on  steps  determined  by  the  least 
work  in  the  breaking  process. 

It  is  clear  from  (60)  that  the  distributions  of  the  Qj  depend  on  the  Uj.  We  know  the  Uj 
depend  on  e,  z,  T,  J0  or  J0  and  As,  and  h(i).  It  follows  that  the  distributions  of  the  Qj  depend 
on  the  same  quantities.  Thus,  it  is  possible  to  study  how  changes  in  these  quantities  influence 
the  distributions  of  the  Qj.  However,  we  will  not  attempt  a  study  of  this  point  in  this  paper.. 

There  are  a  number  of  additional  parameters  that  must  be  estimated  if  the  model  is  to 
account  for  the  random  nature  of  the  Qj  or  Zy  It  appears  that  the  same  high  replication  tests 
needed  to  estimate  the  material  parameters  in  the  deterministic  model  will  suffice  to  supply 
estimates  for  these  additional  parameters.  Thus,  the  same  set  of  tests,  at  the  same  e ,  z,  and  h(x), 
is  all  that  is  needed  to  render  the  model  a  predictive  tool  whether  or  not  randomness  is  included. 

Finally,  we  mention  that  [22]  contains  some  of  our  initial  thoughts  on  how  reaction  rate 
theory  might  be  applied  to  the  FCG  problem. 
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7.  Conclusion 

We  present  a  deterministic  macro  model  of  FCG  based  upon  a  micro  result  provided  by 
reaction  rate  theory.  A  center  crack  panel  under  periodic  tensile  load  is  the  basic  physical 
situation  considered.  The  model  contains  explicitly  the  stress  distribution  s(t,q)  on  the  unbroken 
portion  of  the  crack  plane,  the  for  field  stress  sQ(t),  the  absolute  temperature  T,  and  two  material 
parameters  AG~  and  (3  of  which  AG~  can  be  estimated  from  available  information,  but  the 
growth  control  parameter  [3  must  be  estimated  from  test  data  since  it  depends  on  how  the 
material  is  formed  into  specimens/components  and  the  environment. 

For  a  reasonable  choice  of  s(t,q)  that  assumes  no  dynamical  effects  present  in  the 
specimen/component,  we  find  an  evolutionary  equation  in  the  dimensionless  crack  length  y  and 
two  other  quantities  5  and  ra  that  describe  aspects  of  the  evolution  of  the  stress  distribution  Joint 
evolutionary  equations  for  these  quantities  are  needed  to  determine  y,  5,  and  rc  as  functions  of 
time,  and  to  determine  if  they  form  a  state  vector.  We  conclude  from  this  approach  that  the  FCG 
process  is  more  complex  than  can  be  described  by  y  alone. 

At  the  present  time,  no  such  joint  evolutionary  equations  are  available.  To  proceed  with  a 
study  of  some  of  the  capabilities  of  the  model  based  upon  a  statistical  analysis  of  curiendy 
available  data,  we  introduce  constraint  equations  for  5  and  ra  that  reflect  general  properties  of 
polvcrystaline  ductile  metals  in  fatigue,  and  use  the  resulting  deterministic  constrained  equation 
to  describe  the  evolution  of  v,  containing  the  two  constraint  parameters  e  and  z  in  addition  to 
AG4-,  (3,  T,  and  the  far  field  tensile  load  s0(t)  =  sj,  +  As  h(wt),  where  h(t)  is  the  periodic  wave 
form.  An  equation  of  this  type  only  can  apply  to  the  mean  behavior  of  the  crack  growth. 

Available  reasonably  high  replication  FCG  data  are  for  simple  harmonic  tensile  load 
variation  and  specimens  of  polvcrystaline  ductile  metals.  For  a  load  of  this  type,  the  constrained 
equation  for  the  evolution  of  v  can  be  useful  approximated  by  a  difference  equation  in  discrete 
time  and  discrete  crack  length  that  not  only  makes  parameter  estimation  of  [3  straight  forward, 
but  also  permits  us  to  study  the  effect  of  changes  in  T,  h(t),  ?0,  As,  and  [3  on  the  FCG  process. 
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We  conclude  from  a  parameter  estimation  study  based  upon  these  data  and  the  mean 
model  that 

a)  there  can  be  transient  behavior  in  the  (3j70  estimates  due  to  a  temperature  change  or  a  load 
change, 

b)  the  growth  control  parameter  (3y  depends  on  the  crack  length  through  the  relation 
Py  =  P0(l-cy)  once  the  transient  effect  of  the  initial  temperature  change  ceases, 

c)  there  are  three  material  parameters  AG+,  P0,  and  c  when  there  is  only  one  observable  y  and 
a  fixed  set  of  constraint  parameters  e  and  z, 

d)  the  relation  Py  =  p0(l-cv)  is  due  to  the  work  hardening  properties  of  the  material, 

e)  P0  depends  on  70  and  As  but  c  does  not  appear  to  do  so,  both  depend  on  e  and  z,  p0 
depends  on  T  but  a0  and  c  do  not  appear  to  do  so,  p0  and  c  may  depend  on  h(t),  and 

f)  two  or  three  suitable  tests  are  needed  to  estimate  PD  and  c  so  that  the  model  can  predict 
crack  growth  behavior  over  a  reasonable  range  of  70  and  As  values. 

Further,  we  conclude  from  e)  that  when  s0(t)  is  oscillatory  but  not  periodic,  transient  material 
behavior  may  occur. 

The  constrained  model  or  difference  equation  in  Avj  and  Atj  contains  explicitly  the  wave 
form  h(x),  T,  70  and  As.  We  conclude  that  predicted  Atj  for  a  given  Ayj  is  sensitive  to  changes 
T,  Sq  and  As,  and  the  predicted  Atj  can  change  substantially  for  changes  in  h(r)  with  the  same 
70  ,  As,  and  T.  Thus,  we  conclude  that  when  conducting  experiments  and  when  predicting  crack 
growth  in  service  these  quantities,  particularly  T  and  h(t),  play  a  significant  role. 

The  possibility  of  increasing  only  T  to  accelerate  a  test  is  a  useful  observation  since  such  a 
procedure  appears  to  eliminate  load  interaction  problems  that  may  arise  when  ?0  or  As  are 
changed  with  T  constant,  as  is  usually  done.  In  addition,  we  note  that  this  method  of 
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accelerating  a  test  can  be  used  with  profit  when  conducting  the  two  or  three  tests  needed  to 
determine  how-  po  and  possibly  c  depend  on  the  load  which,  as  we  have  remarked  above,  are 
needed  to  predict  life  under  different  70  and  As  than  employed  in  these  tests.  We  conclude  that 
the  explicit  presence  of  T  in  the  model  is  an  important  feature  of  this  model. 

The  constrained  model  with  one  observable  v  contains  explicitly  the  far  field  stress 
variation  s0(t).  Thus,  by  means  of  numerical  integration,  we  can  predict  crack  growth  behavior 
for  any  tensile  s0(t),  including  static  load.  However,  we  conclude  from  a)  and  e)  above  that 
transient  material  behavior  may  put  accurate  prediction  based  on  this  model  into  question  when 
s0('t)  is  oscillatory  but  not  periodic.  We  conclude  that  for  accurate  prediction  for  any  oscillatory 
non  periodic  -’0(t)  it  may  be  necessary  to  have  the  joint  evolutionary  equations  for  v,  5,  and  ra 
for  it  is  only  by  such  a  model  that  we  may  be  able  to  take  into  account  the  load  history 
dependence  of  the  material. 

The  unit-jump  B-model  [19]  contains  the  probability  transition  matrix  P  with  elements 
pi,  Qj ,  (pj  +  ^  =  1)  on  the  jth  line.  Since  (bj  —  bj_i )/qj  =  E(AT, },  we  conclude  that  it  is  now 
possible  relate  the  (bj  -  by_i  )/qj  to  T,  ?0  ,  As,  and  h(x).  This  is  an  important  point. 

The  model  assumes  that  the  rate  of  single  bond  breaking  Kb  is  determined  by  a  single 
energy  barrier  and  there  is  no  bond  healing.  It  is  clear  that  it  is  possible  to  relax  our  assumptions 
in  order  to  accommodate  other  possibilities.  We  conclude  that  the  simple  model  presented  here 
permits  extension  in  a  number  of  directions  of  this  type.  Bond  healing  can  be  accommodated  in 
B-models  bv  assuming  nonzero  elements  below  the  main  diagonal  in  P. 

The  specimen  geometry  and  loading  can  be  changed.  A  new  model  would  have  to  be 
created  for  each  new  situation.  This  constitutes  another  area  into  which  this  type  of  model  can 
be  extended. 

If  there  are  dynamical  effects  in  the  specimen,  s(t,q)  cannot  have  the  simple  form  we  have 
assumed.  It  may  be  possible  to  include  this  feature  although  we  have  not  considered  it. 

Variability  in  material  behavior  is  introduced  into  the  deterministic  difference  equations 
by  the  simple  expedient  of  multiplying  the  growth  control  parameter  by  (1  +  Qj),  where  the 
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Qj  are  a  family  of  random  variables  indexed  on  (yrl ,  yQ;  Atj  is  then  replaced  by  the  random 

variable  AT,.  We  already  know  from  previous  work  [17,18]  and  also  from  the  statistical  analysis 

conducted  here  that  the  ATj  are  dependent  and  non-normal  with  positive  skewness.  We 

conclude  from  the  analysis  of  the  VUG  and  GD  Data  presented  above  that  the  Qj  are  dependent 

and  non  normal  with  negative  skewness.  We  also  conclude  that  the  distribution  of  the  Qj  and 

their  moments  depend  on  ?0  ,  As,  and  T  for  given  e  and  z,  although  given  the  database  it  is  not 

/ 

possible  at  this  time  to  determine  the  nature  of  this  dependence  on  sQ  and  As.  For  the  same 
reason,  we  can  not  determine  the  nature  of  this  dependence  on  h(t). 

The  closeness  of  the  standardized  edf  s  of  the  Atj  and  of  the  Qj  requires  explanation,  since 
it  suggests  that  there  is  an  underlying  random  material  process  present.  It  is  possible,  as 
mentioned  earlier,  that  if  we  introduce  material  randomness  at  the  micro  level  and  develop  a 
stochastic  macro  model  from  this  starting  point.  We  can  find  an  explanation  for  this  interesting 
experimental  fact. 

Parameter  estimation  for  the  above  stochastic  model  can  be  handled  from  data  needed  to 
estimate  the  parameters  of  the  deterministic  model. 

Based  upon  the  above,  we  conclude  that  the  model  appears  capable  of  illuminating  many 
of  the  complexities  that  are  experimentally  observed  in  the  FCG  process,  it  has  useful  predictive 
possibilities,  and  it  suggests  a  simple  method  for  accelerating  FCG  tests. 
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Figure  Captions 


E24.04.01  Round  Robin  Data.  Estimates  of  pjSo  vs.  crack  length  a. 
a  =  2.55  v,  s  =  1,  z  =  0.25 

GD  Data.  Estimates  of  p;?o  vs.  crack  length  a.  h(x)  =  sin  x,  21  points, 
a  =  50.8  y,  e  =  1,  z  =  0.25. 

/ 

GD  Data.  Estimates  of  po  vs.  ?o.  h(x)  =  sin  x,  e  =  1,  z  =  0.25. 

EHN  Data.  Estimates  of  pjSo  vs.  crack  length  ax  for  variable  e  and  fixed 
z  =  0.25.  h(x)  =  sin  x,  a  =  35  y. 

IHN  Data.  Estimates  of  PjSo  vs.  crack  length  a  for  variable  5  and  fixed 
e  =  1.0.  h(t)  =  sin  x,  a  =  35  y. 

VUG  Data.  Effect  of  changing  T  on  a  vs.  t  curve.  h(x)  =  sin  x,  a  =  76.2  y, 
e  =  1,  z  =  0.25. 

VHG  Data.  Effect  of  changing  h(x)  on  a  vs.  t  curve.  a  =  76.2y,  6=1, 
z  =  0.25. 

VHG  Data.  Smoothed  edfs  of  standardized  qj  for  variable  Aa  and  fixed 
load.  h(t)  =  fin  x,  a  =  76.2  y,  €  =  1,  z  =  0.25. 

VHG  DATA.  Smoothed  edfs  of  standardized  qj  for  fixed  Aa  and  fixed 
load.  h(x)  =  sin  x,  a  =  76.2  y,  e  =  1,  z  =  0.25. 
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Figure  10:  GD  Data.  Edf' s  of  standardized  qj  for  fixed  Aa  and  variable  load. 
h(x)  =  sin  x,  s  =  1,  z  =  0.25,  16  —  17  mm. 

Figure  11:  GD  Data.  Edf’s  of  standardized  qj  for  fixed  Aa  and  variable  load. 
h(x)  =  sin  t,  e  =  1,  z  =  0.25,  17  -  18  mm. 

Figure  12:  GD  Data.  Edf  s  of  standardized  qj  for  fixed  Aa  and  variable  load. 
h(t)  =  sin  x,  e  =  1 ,  z  =  0.25,  16  -  18mm. 


•  • 


0.1  — r 


edf 


-  25-27  (mm) 

-  25.6-26.6  (mm) 

-  25.8-26.2  (mm) 

-  24.6-27.6  (mm) 


(qj-m)  /a 


-3  0  3 


(qj-m)  /c 


i  -=>  * 


o 


